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ABSTRACT 


The system of 13 Ceti, in which three components are now known with certainty, 
has always presented unusual difficulties in observation and interpretation. Several 
orbits, visual as well as spectroscopic, have been derived—none as yet entirely satis- 
factory—and many suggestions have been made concerning various possible sources 
of perturbation. 

It is shown in the present paper that the orbit derived by Maggini from a large 
number of interferometer measures does not obey the law of areas, and from a rediscus- 
sion of the entire material a new visual orbit is derived which accords closely with 
Aitken’s original one. The various suggestions and explanations offered for the observed 
perturbations are discussed, and it is shown that most of these are quantitatively un- 
reasonable, while furthermore it appears that in most cases the perturbations were a 
direct result of the use of too many uncertain decimals. 

The spectroscopic orbit is derived anew from all the available material and it is found 
to possess a definite, though small, eccentricity, while furthermore a large advance of 
periastron appears to be indicated. This is tentatively explained as due to the pre- 
cession operating in the spectroscopic system owing to the large tidal deformation of 
the two components. 

Finally, by making plausible assumptions concerning the orbital planes, values are 
derived for the absolute dimensions of the spectroscopic system and for the individual 
masses of the three components now known. 


1. The star 13 Ceti (a 030™,6—4°09',1900.0,5.6,6.4,Go) was 
discovered to be a close visual double by Hough in 1887. The 
brighter component was found to be a spectroscopic binary by Frost’ 
in 1907, while orbits were first calculated by Aitken’ (visual) in 
1907, and by Fox (spectroscopic)’ in 1908. Paraskevopoulos‘ and 

* Astrophysical Journal, 25, 60, 1907. 

2 Publications of the Lick Observatory, 12, 5-6, 1912. 


3 Astrophysical Journal, 27, 375, 1908. 4 Ibid., 52, 110, 1920. 
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Pogo’ further investigated the spectroscopic system, and Maggini® 
derived a new orbit for the visual pair from a long series of inter- 
ferometer measures. The system is one of great complexity and 
presents unusual difficulties in the observations, especially in the 
spectroscopic observations where only one spectrum, that of the 
brighter component of the spectroscopic pair, appears with certain- 
ty, though it may be disturbed by being blended with that of the 
fainter visual component. In an effort to explain the many discrep- 
ancies present in the observations, numerous suggestions have been 
made. In discussing the orbit of the spectroscopic pair, Paraske- 
vopoulos claimed to have found clear evidence of perturbations in 
it by the fainter, more massive, visual component during periastron 
passage. The material used consisted of eight spectrograms taken 
during the opposition of 1906, ten taken in 1907, eleven in 1908, and 
ten in 1912. The first two sets were combined into one group, and 
for each of the three groups a separate set of elements was derived. 
The ten observations in the last set were subjected to the method of 
least squares and “definitive” elements were calculated. 

Paraskevopoulos then claimed that a comparison of these three 
sets of elements showed clearly that the spectroscopic period is longer 
during visual periastron than during apastron, that the value of a 
sin 7 varies, owing to changes in either a or 7 or both, and that the 
motion of the line of apsides is well established. The three systemic 
velocities which hold for different positions in the visual orbit were 
then used to calculate the true systemic velocity y, = 7.134 km/sec., 
and a constant, which involves the parallax and the visual mass ratio. 
Finally, four more observations made in 1919 were adduced and 
compared with the velocities calculated from the 1912 elements; the 
resulting differences O—C, which run as high as 30 km/sec., were 
taken to show “very clearly the variation of the elements of the 
orbit.” 

Pogo, from much more extensive spectroscopic material, claims 
that the spectroscopic orbit is circular, and attributes such irregulari- 
ties as are found to exist in the observations to changes in the at- 
mospheres of the stars or to disturbing effects arising from the spec- 
troscopic duplicity of the fainter visual component. Maggini sug- 

5 Ibid., 68, 116, 1928. 6 Astronomische Nachrichten, 233, 97, 1928. 
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gests that the fainter visual component B has an invisible companion 
and moves in an orbit with a radius of about 0703-0704. Finally, 
Slavenas,’ in the course of a photometric investigation, used 13 Ceti 
as a comparison star and suspects it to be variable by o™2 with a 
period of 044692. 

2. In making these various suggestions and explanations the fun- 
damental point whether they are quantitatively possible or even 
reasonable seems to have been continuously overlooked. 

In the discussion of Paraskevopoulos, thirty-nine observations in 
all were used to calculate no less than eighteen different unknowns; 
furthermore, the values O-C tabulated for each observation indi- 
cate a mean error in each velocity of about 2.6 km/sec., hence it is 
evident that not more than three, and certainly not seven places are 
warranted in a sin 7, that the differences in K, a sin 7, and e need not 
be real, and that w must be exceedingly uncertain. Where the inter- 
vals covered by the 1908 and 1912 observations measure only some 
sixty and thirty-six periods, respectively, the uncertainty in the pe- 
riod may well be as much as o4o1 and hence even the differences in 
period need not be real. With as smal] an eccentricity as the present, 
it is certainly not warranted to give T to o40001, especially since no 
mention has been made as to whether or not the observations haye 
been reduced to the sun. Of the eleven observations in 1908 no less 
than nine deviate systematically from the curve drawn, while the 
large values for O—C for the four isolated observations of 1919, in- 
stead of ‘‘showing clearly” that perturbations exist, would, in the 
writer's opinion, furnish additional evidence that the calculated ele- 
ments are uncertain. The least-squares solution which gave yo+ 
7.134 would, if its precepts had been followed, yield for the mean 
error of y, + 2.3, and for C, 11.372+4.3. 

The evidence presented by Paraskevopoulos does not, therefore, 
prove that perturbations exist, but indicates rather clearly that the 
material is insufficient to draw definite conclusions. 

It is true that, qualitatively at least, the perturbation of the period 
is in the right direction. Quantitatively we may, perhaps, make a 
very rough estimate by borrowing some results from the lunar the- 
ory. Asa first approximation we may neglect the eccentricity of the 


7 Ibid., p. 125. 
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spectroscopic orbit, and assume it and the visual orbit to be co- 
planar. There then is no orthogonal perturbation, and the chief 
perturbative force active in lengthening the period is the normal 
component whose average value is 
e 
— 2S. p, $ 
where, in the lunar theory & is the Gauss constant, S the mass of the 
sun, r the distance earth-moon, and R& the distance earth-sun. The 
ratio € of this to the attractive force of the earth on the moon is 


S r3 


2(E+M) Rs’ 


the consequent lengthening of the period, since 1/P © 1 mass, would 
be roughly given by 
nr. ss 3 
Pp * 4(E+M) Rs" 

In the case of the sun-moon, S/(E+M) is 3.2105 and R/r=4oo, 
hence AP/P =1/800, which is at least of the right order of magni- 
tude, as the actual lengthening is about one hour, or 1/650.° 

As we shall see later, we may reasonably assume in the present 
case Mg=I.1, M4,44,=1.3, F=0.035 astronomical units, R=1.29 
A.U. at visual periastron, and 8.5 A.U. at visual apastron. We then 
find AP=o0.8 X10 5 days at periastron, and, of course, a negligible 
quantity at apastron. Paraskevopoulos finds AP =0.4 X 10 *, or five 
hundred times as much, and while undoubtedly the simple lunar 
theory is not applicable in this instance, where the mass ratio A,/ A, 
is about 3:1 instead of 81:1, it is hardly likely that this would in- 
crease the effect five hundred fold. A further, more serious objection 
is that Paraskevopoulos has overlooked the fact that his longest 
spectroscopic period corresponds to the mean epoch 1912.9, or o.1P, 
after periastron passage. In an orbit as highly eccentric as the pres- 
ent visual one, this gives R=4.1 A.U., or nearly half the apastron 
distance, and hence AP only 0.3 X10 ° days, or 1/15,000 of that 


8 Moulton, Celestial Mechanics, p. 349, 1928. 9 Tbid., p. 348. 
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claimed by Paraskevopoulos. We are thus forced to the conclusion 
that the actual perturbations are much too small to be observable. 

Pogo’s claim that the spectroscopic orbit is circular will be taken 
up farther on, and shown to be erroneous. 

Following up Maggini’s suggestion, we find that the semiaxis 
major of the relative orbit of B, and B, is at least o%07, or equal to 
the periastron distance of A and B, hence one would expect severe 
perturbations to result. 

Furthermore, it is seen that although Maggini’s ninty-four ob- 
servations are extremely concordant, and define the apparent ellipse 

TABLE I 


MAGGINI’S INTERFEROMETER MEASURES 


Dat P Area Per Cent At Areal 

é 4 r . euene , . 

une of Ellipse Years Constant 
» °o ” } 

1922.95 283.5 0.200, | me 

1924.08 306.0 .208} | $ 7 — = = 

1925.19 339.8 | .126! ; 

1925.62 20. .060 | : 

1926.00. 171. 070} 7” = -083 

1926.10 19O. .1054 | 

1920.20 198.6 | .130> 166 | 0.82 . 202 

1920.92 228.7 . 2304 aie 0.03 | 178 

1927.85 249.4 278} coe ee | . sue 

1929.84.... 283.5 0.260! baits “99 Weee 


with great precision, while also this apparent ellipse fits the visual 
observations tolerably well, the law of areas is not obeyed in this 
ellipse by the interferometer measures themselves. Taking several 
normal places whose positions appear to be indicated with almost 
absolute precision from the observations, we find the values shown 
in Table I for the areas described (expressed as a fraction of the 
whole ellipse) and for the times consumed. 

It is evident that the apparent periastron distance is too small, 
hence the eccentricity is too large—and this is the chief discrepancy 
between Maggini’s orbit and Aitken’s visual one. Also, it must be 
regrettably admitted that the interferometer measures are apparent- 
ly subject to large systematic errors. (If Maggini’s orbit is aban- 
doned, his suggestion concerning the duplicity of B automatically 
follows. ) 
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Slavenas’ suggestion of a variation in total light of o™2 with a 
period of 044692 is very hard to fit in with any data at the present 
time. The length of period suggests either a cluster variable or an 
eclipsing binary. Since the brightest star of the system is a G dwarf. 
the first suggestion is ruled out. The spectroscopic period of A is 
2.08, hence it would seem that only B could be the eclipsing binary, 
but if the total light of AB changes by o.2, then the fluctuation of B 
alone should surpass o™5. It should be possible to verify this by vis- 


TABLE II 


ELEMENTS OF THE VISUAL ORBIT OF 13 CETI 











VISUAL INTERFEROMETER 

Aitken | Present Paper | Maggini acon 
| | Corrected 

7 Re: <2) ee 28°7 2°5 23°0 
1* +53.4 | +53.0 +52.6 | +50.0 
w 66.8.. | 75.0 | 85.0 | 65.0 
a 0.242 | 0.249 | 0.247 | ©. 24 
€ 0.725 Were 0.730 | 0.816 | 0.71 
P Rr ok, eee | 6.93 | 6.895 | 6.92 
(OR ee yy re I9:9.20 | 1925.92 | 1925.04 








* The positive sign for 7 was assumed to be known from spectroscopic observations. 


FINAL ELEMENTS 


Q= 24° ¢= 0.73 
1=+52 r= 6.91 years 
w= 73 T= 1925.94 
a=o"24 Angles increasing 


ual estimates, but until it is done the eclipsing binary character of B 
cannot be considered as established. 

Finally, the visual orbit derived by Aitken shows residuals O—C 
with the observations of as much as 4o’ in angle. 

In view of the generally unsatisfactory state of the spectroscopic 
as well as the visual system, it has been thought worth while to dis- 
cuss all available observations anew. In what follows the visual 
components will be referred to as A and B, the spectroscopic com- 
ponents of A as A, and A,. 
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3. Visual orbit.—Using all available visual observations given in 
Burnham’s and Aitken’s catalogues, also a recent observation by 
Van den Bos, new elements were calculated. In this connection I 
am much indebted to Mr. Arnold Rosenstein, a graduate student in 
our department, who performed all the actual calculations. It 
should be emphasized that his calculations were made entirely with- 
out knowledge of Aitken’s or Maggini’s orbit, and hence the results 
should be wholly independent. 140 
In Table II are given, in order, 
Aitken’s elements, those derived 
here, Maggini’s interferometer 
orbit, and a set of elements de- 
rived from Maggini’s observa- 
tions by roughly reshaping the 
ellipse to make it fit the law of 
areas. The values finally adopted 
for the elements are given below 





x Visvac 


O /wreepeeome zee 


the table; the graphical repre- O 
sentation of the observations is Fic. 1.—The visual orbit of 13 Ceti. 
shown in Figure I. Crosses represent visual measures; open 


circles, normal places derived from Mag- 


Two trigonometric values for 
8 gini’s interferometer measures (cf. Table I). 


the parallax are available; viz., 

0”048 (S;, rel.) and o%051 (Y,,; rel). Using the systematic correc- 
tions derived by the writer’ and taking into account the spectro- 
scopic parallax, we would obtain o’o50 for the absolute value of p. 


Hence the absolute dimensions of the system are 


a=4.9 astronomical units, g=1.3 astronomical units, 
total mass=2.4 ©. 


4. Spectroscopic orbit—Using the new visual elements, we may 
now calculate the relative radial velocity of the two components, the 
systemic velocity of A and the corrections to be applied for light- 
time at various epochs due to the orbital motion of A around the 
center of gravity of AB. Actually observed systemic velocities of A 
are available at many different dates, and may be taken directly 


10 Proceedings of the National Academy of Sciences, 16, 472, 1930. 
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from Pogo’s paper. At this point it is a pleasure to acknowledge the 
writer’s obligation to Dr. Otto Struve for calling attention to this 
paper, of which the writer had been totally unaware. The much 
larger spectroscopic material that has become available through 
Pogo’s paper has made it possible to verify the conclusions derived 
to a much greater extent than before. 

If I’ be the true systemic velocity of AB, V the orbital velocity of 
B with respect to A, y the observed systemic velocity of A, and £ the 
ratio: mass of B/mass of total system, we have 


y=T—-B-V. 


Performing the calculations for those epochs for which Pogo has 
given the observed systemic velocities of A, and which rest on five or 
more observations,'' we obtain the data collected in Table ILI; these 
are plotted in Figure 2. 


TABLE Ul 
t V , t V ¥ 
1g06.9go. =I =FIOVO TOELT3 O05 5 « 6 =FI7.0 
07.87 — 6 -FII.O 22 :592....505 A aeESs4 
08 . 89 —1z + 8.0 25.95..+..0.$29 + 7.0 
EEs72 —=25 — 320 BOOT Sos os —16 +13.0 
12.97 —216 +17.0 47 BLS | gO se = 


The peculiar differences, large in amount and systematic in char- 
acter, already mentioned by Pogo, are clearly seen. The slope of the 
points, which defines 8, however, is unmistakable, and we may per- 
haps adopt for the present the graphically derived values 


[=+11.0; B=0.45 


and hence the following values for the masses A: 1.27, B:1.05, the 
second decimal being added merely for purposes of calculation. 

Using 8, it is further possible to calculate the corrections for light- 
time due to the orbital motion of A; making also the usual heliocen- 
tric light-time correction, we may now reduce all spectroscopic ob- 
servations to a uniform system. 


"Op. cit., Table XI. 
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All observations given by Paraskevopoulos were first dealt with; 
the four sets made during the oppositions of 1906, 1907, 1908, and 
1912 were reduced separately with a preliminary value of 24082 for 
the period, the times of maximum and minimum velocity read off, 
and these used to calculate a second approximation for the period, 
viz., 2408186, with a mean 
error of about three in the 
last place. With this all ° 
observations were reduced . 
to the mean epoch of all, 
around J.D. 2418402; fur- 7 
thermore, the velocities 
were corrected for the cal- 
culated systemic velocity of Fic. 2.—Abscissae: relative velocity of B 
A—that due to its orbital  ‘"eferred to A, as calculated from the visual 
A aa orbit; ordinates: observed systemic velocities of 
motion, as well as the ‘‘true Lito 
systemic velocity of AB. 

The plot of all thirty-nine observations is given in Figure 3. Con- 
sidering the difficulties and uncertainties inherent in the observa- 
tions, the velocity-curve is well defined; the four isolated observa- 
tions made in 1919 have been added, but have not otherwise been 
used. Planimeter readings give a residual systemic velocity of +2 
km/sec., which was to be expected considering the distribution of 
points in Figure 2. 

The outstanding feature of the curve is the great sharpness of 
minimum as contrasted against the relative flatness at maximum. 
Coupled with the fact that the curve is practically symmetrical 
around the ordinates of maximum and minimum velocity, this shows 
clearly that the orbit is eccentric, and that the longitude of perias- 
tron is roughly 180°; T is therefore immediately given as the instant 





i 


4420 K m / See 











V 0 +20 Kw/sec 


of minimum velocity. 

The same procedure cannot be immediately applied to Pogo’s 
forty-seven observations made during 1923-1928 since it is evident 
here that the systematic correction to systemic velocity is not the 
same in each year. Since for the present we are concerned only with 
the shape of the velocity-curve, it is perhaps allowable to determine 
these corrections empirically. A reasonably good agreement may 
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be obtained by applying the following corrections to the veloci- 
ties. 


Ce a ...—9 km/sec. 
NOQUO isos, Soe hk BR oe ee EO ; 
Cob feck (0 |: i a es 


The resultant velocities were then reduced to the mean epoch of 
all, around J.D. 2424548. Again it is seen that they conform toler- 
ably well to a smooth curve, which is plotted in Figure 4. 





+20 


-40 














1 
2.00 JSD. 2,4/8,403.00 


Fic. 3.—The velocity-curve of 13 Ceti A during the years 1906-1912. Abscissae 
indicate the phase, reduced to the common epoch of J.D. 2418402, after correction for 
heliocentric and astrocentric light-time had been made. Ordinates indicate the observed 
velocities, corrected for true systemic velocity of AB and orbital motion of A. 


Calculating the elements of the orbits corresponding to Figures 3 
and 4 in the usual way, we obtain the data given in Table IV; the 
system velocities have, of course, been omitted since these now meas- 
ure only the amount of systematic error in the observations. 

The close agreement in the value of e is doubtless accidental, but 
in any case it is clear that the small values of 0.015 + 0.018 and 0.022 
+o0.018 obtained by Pogo for the 1924 and 1928 curves are entirely 
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incompatible with the marked eccentricity of the 1906-1912 curve. 
Pogo further remarks on the change in w found by Paraskevopoulos 
that not only do we “‘have the case of a circular orbit, not the case of 


















+4, 
+40 4 “= 4 
+20 + 4 
P 
°o — al 
+ /923 
- /924 
A 1925 
~oa.> @ /926 4 
© /927/28 
-4Ao} (oY) a a 
l l 
J.D. 2,424,548.00 Goo 


Fic. 4.—The velocity-curve of 13 Ceti A during the years 1923-1928. Abscissae 
indicate the phase, reduced to the common epoch J.D. 2424548, after correction for 
heliocentric and astrocentric light-time had been made. Ordinates indicate the observed 
velocities corrected for orbital motion of A, systemic velocity of AB, and such empirical 
corrections as appeared necessary to reduce them to the same level. 


the motion of the line of apsides’” but also that “there is another 
reason why it would be futile to interpret these discordant results as 


TABLE IV 


1906-1912 1923-1928 
|: GPR aries eye 36.8 38.4 
Oi Se ce We 0.10 O.II 
@. | Sexo eee 180. 280. 
f (pF Re, Ae J.D. 2418402. 33 J.D. 2424548 .60 
asimi.....:.. ‘.egdtse 1.08 X 10° 


P= 2'08186 


a motion of the line of apsides.’”’ The reason advanced is that the 
perturbation produced by the fainter visual companion would 
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change the direction of the line of apsides by something like 1° in 
several hundred years, according to a formula derived by Slavenas.” 

While the present writer agrees with this statement in so far as it 
applies to the orbits derived by Paraskevopoulos and the perturba- 
tion caused by B, the orbit can most certainly not be dismissed as 
circular at the present time. Furthermore, in the present writer’s 
opinion the systematic uncertainties in the systemic velocities make 
it futile to apply a least-squares solution in order to calculate the 
small but meaningless probable errors in e. Finally, blaming any 
possible motion of the line of apsides on the companion B is looking in 
the wrong place for there is another far more likely cause for such 
motion. 

Using the masses for A and B determined above, viz., 1.27 and 
1.05, respectively, we find, with P= 2408 =o0%0057, a (spectr.) = 
0.035 astronomical units. Since a sin 7 for the absolute orbit of A, 
is 1.06X10° km, or 0.0071 A.U., the mass ratio A,/A, cannot be 
more than 4:1. Making the reasonable assumption that the visual 
and spectroscopic systems are approximately coplanar, we have sin 
i1=0.8 and a (A,) =0.009 A.U., giving a mass ratio of 3:1 in the spec- 
troscopic system and individual masses of 0.95 and 0.32, respective- 
ly. Since the spectral class of A, is G and its absolute magnitude 4.1, 
it would thus seem that this brighter component is very similar to 
the sun, while A, might well resemble an M dwarf; hence their diam- 
eters should be of the order of 1.5 X 10° and 0.7 X 10° km, respective- 
ly, i.e., the separation of their surfaces is only about three to four 
times the sum of their radii. Presumably the axial rotations of both 
components are synchronized with their revolution; this would make 
the ratio ¢ of the centrifugal force at the equator of the larger star to 
the attractive force about the same as for the earth. But where the 
two stars are so close their tidal deformations may well be large, and 
hence we should theoretically expect a rapid advance of the line of 
apsides. 

Even though all numerical quantities involved are extremely un- 
certain, it might be of value to see whether the observed motion of 
the line of apsides appears reasonable. Where the fainter component 


21 bid., p..133- 
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of the spectroscopic system is presumably of less mass, of much 
smaller size, and of considerably higher density than the primary, we 
may perhaps, as a first approximation, treat it as a satellite revolving 
around an oblate planet. In that case the ratio K of the period of ro- 
tation of the line of apsides to the period of orbital revolution is 
given by 


where a is the semiaxis major of the orbit, 7 the radius of the larger 
stars, € its oblateness, and ¢ the ratio mentioned above. In the pres- 
ent case we have K = 104,a=5.2X10° km, r=7 X10 km, g= 1/300; 
hence €= 1/150, which appears not unreasonable. 

A glance at the data of Table IV shows that periastron has ad- 
vanced by about 1oo° during the 2952 periods which have elapsed 
between the two sets of elements, but this number of 100 may well 
be uncertain by some 20°-30°. Hence we find for the rate of advance 
of periastron 0.03°/P. Since a true anomaly of 100° corresponds, for 
e=o.10 and P= 2408, to a time interval of of51, we finally obtain 


Anomalistic period. ............. . 2408207 
Siebetenh MATION. oss ceo beechans 2408189 


With the adopted velocity of +11 km/sec. for the center of mass 
of the entire system, this gives for the “true” sidereal period, freed 
from Doppler effect, P = 2408182. 

In order to settle the matter conclusively, it would probably be 
necessary to obtain a large number of observations during at least 
two more epochs separated by some ten years. At the next visual 
apastron passage it might perhaps be possible, by setting the slit of 
the spectroscope at right angles to the visual position angle, and by 
using only nights of excellent definition, to obtain some impression 
of the spectrum of the fainter visual component, and to decide 
whether this component, too, is perhaps not a spectroscopic binary, 
which seems not improbable on account of its mass. 
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5. Summary.—Tentatively combining the two spectroscopic or- 
bits into one, we may give the following summing-up of our present 
knowledge concerning the system 13 Ceti: 


R.A. 0:30 Dec.—4:09 5.6 Go, 6.4—p=0"o050 
w=o"41 in 93°, or T=39 km/sec., V=+11 km/sec. 


Visual Orbit Spectroscopic Orbit 
P=6.91 years P= 2.08189 days 
ad=0"244=4.9 astr. units a sin i=1.06X10° km 
e€=0.73 e=0.1 
Q= 24° K=37 km/sec. 
i=+52° 
w= 73° w=180f0°o16(t—T) 
T= 1925.94 T=J.D. 2418402. 33 


Angles increasing 


MassEs: A,=0.95, A.=0.32, B=1.05 


MINNEAPOLIS, MINN. 
June 14, 1933 


OPACITY FORMULAE AND STELLAR LINE 
INTENSITIES' 
By HENRY NORRIS RUSSELL? 


ABSTRACT 


Theory.—Formulae for the general opacity of a stellar atmosphere are reviewed, 
especially Pannekoek’s approximations. General equations are derived (§ 5) for the 
numbers of neutral and ionized atoms and the optical depth in an atmosphere contain- 
ing two partially ionized elements and others not ionized, and simplified by the introduc- 
tion of functions f, which are tabulated (§ 8). Solutions in which only one component is 
supposed to be partially ionized—the rest being completely ionized or completely neutral— 
are sufficient in almost all practical cases. The influence of variations in abundance, 
gravity, and temperature are discussed ($§ 9-10). If the opacity depends wholly on the 
interaction of electrons and ions, the lines of any element will have a maximum (for 
varying temperature) only if some other element of easier ionization is present, and the 
position of the maximum will depend on the relative abundance of the two (§§ 13-16). 
Opacity due to interaction of neutral atoms and electrons produces similar effects (§ 17). 
Equations are given for both cases and for the combination of the two. 

A pplications.—With the aid of empirical expressions for the surface gravity in stars 
of different temperatures (main sequence and giants) [§ 18], the conditions for maxima 
of lines of H, Mgt, and Fe* are worked out for atmospheres composed mainly of hydrogen 
($§ 19-20) and of easily ionized atoms like sodium (§ 21). The observed low temperature 
for Fe* (class F5) and the persistence of arc lines of Fe in the hotter stars can be explained 
in the former case, but not in the latter. Numerical values for the pressure are derived 
from the writer’s data for the solar atmosphere (§ 22). From the maxima of Fet and 
Ti* it follows (§ 23) that hydrogen is from 1000 to 2000 times more abundant than all the 
metals together. Fora model atmosphere, containing 1 atom of K to 3 of Na, 12 of Fe, and 
16,000 of H, the electron pressure and amount of matter above the photos phere are calculated 
for main-sequence stars and giants (Figs. i, 2; § 25). The numbers of atoms engaged in 
producing given spectral lines are then computed, approximate account being taken of 
the transition probabilities (§ 26; Figs. 3-8). The computed temperatures of maximum 
agree closely with observation, including the previously anomalous (K) line of Cat. A 
moderate disagreement for H and Mg* may be explained by a contribution of neutral 
atoms to opacity (Figs. 6-7), but more probably arises from distortion of the energy- 
curve in stars of class A, making their color temperatures too high. The outstanding dis- 
cordance in the absolute-magnitude effect for H and He at high temperatures is attrib- 
uted to the Stark effect. The still more anomalous strength of the Balmer series and 
the enhanced lines in red giants is ascribed tentatively to large differences in tempera- 
ture within the atmosphere. 


1. It is generally recognized that the intensity of lines in stellar 
spectra depends on the depth to which we can ‘‘see down” into the 
outer layers, and hence on the law of opacity. 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 477. 


2 Research Associate of the Mount Wilson Observatory, Carnegie Institution of 
Washington. 
3 E. A. Milne, Monthly Notices of the Royal Astronomical Society, 89, 17 and 157, 1928; 
W. H. McCrea, ibid., 91, 836, 1931; E. A. Milne and S. Chandrasekhar, ibid., 92, 150, 
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The principal source of general opacity is the interaction of elec- 
trons and ions.‘ Part of this arises from ‘“‘free-free’”’ transitions 
which, according to Pannekoek, produce an atomic absorption co- 
efficient for radiation of frequency v 

PFaZ? 


ao(v) =C, 2 ’ (1) 
al. 


where 


6 


_ 4) On ¢ 


g = 2.053 X10 , (2) 
9 he(mk)* 


i 


m being the mass of the electron, P the electron pressure, Z the 
charge of the ion, a the fraction of the total mass which is formed by 
ions of this sort, and F a numerical factor. The equation has been 
derived only for “hydrogenic’”’ atoms, but is probably a tolerable 
approximation for others—especially for the relatively small energy 
changes involved in absorption in the visible region. 

The “‘bound-free”’ absorption due to photoelectric ionization of a 
given atom or ion increases abruptly just beyond the limit of each 
series in its spectrum. Smoothing these irregularities by an approxi- 
mate method, Pannekoek derives an equation which may be written 

PFs 2 


a:(vy)=C, — 
T*y3 hy, 


w 
~~ 


\(hv—kT et *TLERT) , ( 


where F’ is a new constant, jy, the energy of ionization, Z the 
charge, and a the abundance of atoms in which the ionization has 
occurred. Since the atom has here also been assumed to be “hydro- 


1932. Some results of the present investigation were reported to the American Astro- 
nomical Society on December 30, 1931 (see Publications of the American Astronomica! 
Society, 7, 99, 1932). 

4 J. A. Gaunt, Philosophical Transactions of the Royal Society of London, A, 229, 163, 
1930; A. Pannekoek, Monthly Notices of the Royal Astronomical Society, 91, 162, 1931; 
S. Chandrasekhar, Proceedings of the Royal Society of London, A, 135, 472, 1932; B. 
Strémgren, Zeitschrift fiir Astrophysik, 4, 119, 1932; A. S. Eddington, Monthly Notices of 
the Royal Astronomical Society, 92, 364, 1932; R. C. Majumdar, Astronomische Nach- 


richten, 247, 217, 1932. 
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genic,” we may set /v,=Z’x., where x, is the energy of ionization of 
hydrogen. 

To this approximation atoms of all degrees of ionization produce 
the same bound-free absorption. The reason is simple. Absorption 
of this type is produced only by atoms so highly excited that the 
energy x, required to complete their ionization is less than hy. The 
relative numbers of such atoms and of those in the next stage of 
ionization are determined by equations of the usual type, in which 
Xr appears instead of the principal ionization potential. They are, 
therefore, independent of the latter. This formula is inapplicable be- 
yond the principal series limit; but all wave-lengths of astrophysical 
importance are on the safe side. For a mixture, in comparable pro- 
portions, of many elements with overlapping limits, the approxima- 
tion should be good. For hydrogen, Pannekoek’s equations® give 

PF'a x 2 ai 
aslep eee Tips Pac yi ee (4) 
Between \ 3646 and \ 8203 we must set m=3. The ratio, S, of the 
sum to its first term at various temperatures is 
T = 2900° 5800° 8700° 11,600° 23,200 Hn 
S=1.018 1.078 1.128 1.165 1.239 1.368 
The abrupt changes at the series limits have been discussed by 
McCrea.° 
2. It is convenient to express the energies in electron-volts, setting? 





300Xo0 
I,= =13.54, 

€ 
, 300hv 1.234X104% | 
y = 500m = 254X10" | (3) 
: kT eset 
ee a Fe 

e 11,000 ) 


If A is the atomic weight and M, the mass corresponding to unit 
atomic weight,’ we have for the mass coefficient of absorption 

5 Op. cit., p. 163 (bottom of page). 6 Thid., p. 850, 1931. 

7 There ought to be a standard notation for the useful constant here called v. 

8 Birge’s notation (Physical Review Supplement, 1, 62, 1929). Most writers call this 
M jz, but it is not the mass of the hydrogen atom. 
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K(r)=a(v)/AM,. Also, if € is the fraction by number of all the 
atoms present which belong to the given element, x, the fraction of 
these which are singly ionized, and yw the mean atomic weight (ex- 
cluding the free electrons) in grams, au=ex,AM,. Hf, then, we set 


a] Ge 
—_— —_— — 
C1= ay, = 1-27 X 10 , (6) 


we find easily for free electrons and singly ionized atoms 


. 3 
Ko(v) =, cin sis FI, ’ (7) 
a 


and for ionization of neutral atoms 


- 3 
iwc roe. 2. (8) 
a 
where 
G’(A, T)=(V—vT )e" "T+ 0T . (9) 


The factors F’ are well determined theoretically and, under ordinary 
conditions, lie between o.9 and 1.0. We may set F’ =1 without seri- 
ouserror. The F’’sare larger and have not yet been accurately deter- 
mined. Since, according to Y. Sugiura’, * = 6 approximately, we may 
set FI, =80. 

The contributions of different atoms to the opacity are additive. 
Summing for all elements and degrees of ionization, we find 


x(v) =Cyw' PMT "1G", T) ied ox +FI. > eorex) , (10) 


where the sums extend from r=1 upward. The fractions of all the 
atoms which are neutral, singly ionized, etc., are 


Lo= ) €Xo B= Dien... ., 
so that we have 
Kv) =Cu- MPT NGA, T) Det FL Pe} - 


9 Scientific Papers of the Institute of Physical and Chemical Research, Tokyo, 17, 80, 
1931. 


— 


—— ee 


EE 
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When second ionization is negligible, we may write 


«(v) =Cyw MPT #3, G(d, T) ’ (11) 


where 
Ga, T)=(V—oT )e"*T + oT +FI, . 


Hydrogen is an exception. From (4) we find for the visible region 


k(v) =Cyy )3PT “ex, H(T) ’ (12) 
where 
60S = 
H(T) =" 00° eso Fly (13) 


Typical values of G and H are given in Table I. 




















TABLE I 
| 
| G(A,T) 
| Zee 
1 } or | a) | | | | 

| \v= 2.0 | 2.5 3.0 | 3.5 

| A= 6168 | 4034 4112 | 3515 
Meili cetak cin scale st | 290°} 2710 | 5310 | 4.9X104 | 4.5X105 | 3.9X108 
0.50. . | 5800 | 150 | 163 | 377. | +1088 | 3370. 
O2755 . 8700 98 | 99 (| 130 | 204 | 373 
1.00.. 11600 | 88 | 88 99 | 121 164 
2.00.. 23 200 | 82 | 82 84 87 | gI 


3. The mean absorption for the whole flux of radiation follows a 
quite different law. The mean wave-length varies as 7~'; hence the 
integrated absorption varies as PT? for the free-free transitions and 
nearly as PT for the bound-free. The Rosseland mean, which fol- 
lows this law, is to be used in dealing with the general radiative 
equilibrium, but when discussing the optical depth for a given region 
of the spectrum, the monochromatic absorption should be employed. 

The writer is unable to agree in this respect with Milne and 
Chandrasekhar, who use the mean absorption. The very different de- 
pendence on temperature produces important differences in the re- 
sults, as will be seen below. 

At low temperatures, when the ions become few, encounters be- 
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tween free electrons and neutral atoms may have to be considered. 
We may assume, by analogy, 


k2(v) =C3PT~ * yz, . (14) 


Whether C, is even roughly constant as the temperature and com- 
position vary is unknown; but it is probably small in comparison 
with C,. In the cases where this type of opacity is sensible, we may 
set % =1, which is equivalent to adding a small part of the opacity 
due to the ions. 

When ionization ceases at low temperatures, all these forms of 
opacity disappear; but the gas will not actually be perfectly trans- 
parent. Rayleigh scattering at least will remain and give a very 
small value of x, independent of the density. 

Continuous absorption of the type connected with band spectra 
may be more important at the lowest stellar temperatures. This 
should give x proportional to the gas pressure, when compounds are 
just forming. In the absence of detailed knowledge, we may take the 
simple relation 

K3(v) =Cyu , (15) 


where C; may vary with A, and from atom to atom. 

4. We have now to find the strength of spectral lines according to 
these opacity laws. Milne has shown” that a surprisingly good ap- 
proximation may be obtained by treating a stellar atmosphere as if 
it were an isothermal layer of gas, above an opaque photosphere lo- 
cated at a level corresponding to the optical depth 7, for radiations 
near the spectral lines but outside them. The actual lines will be 
nearly equal to those which would be produced by an atmosphere 
transparent except for selective absorption, and of this depth. 

In such an atmosphere the numbers of neutral and singly ionized 
atoms above the level at which the total pressure is p are given by 
the equations 


b 
weNomef Xodp , (16) 


p 
naif xidp . (17) 


© Op. cil., p. 3. The greater the value of xv/x, the better this approximation. 
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Summing for all stages of ionization, we have 
ugn =ep (18) 


for the whole number of atoms of the given kind. The optical depth 
to this level is given by 


iin { "eloddp (x9) 


The equation of ionization may be written 


Sok Syl , (20) 
. : - 
K.= 5 The-H,hT tax) 


The partition functions B are defined by 
B= yi ge-= oT 


where g is the statistical weight and E the excitation potential, and 
the summation extends over all the energy states of the neutral (or 
the ionized) atom. The ratio B,/B, often differs considerably from 
that of the weights of the lowest states." 

In thermodynamic equilibrium 


_2(2mm)'k? _ 


Pm 0.666 . 


C 
Pannekoek has shown” that in an atmosphere of temperature 7” 
exposed to radiation from a photosphere of effective temperature 7, 
equations (20) and (21) may still be used if the photospheric tem- 
perature is introduced and C is multiplied by 3(7’/T)?. Adopting 
T’ =0.86T, after Eddington," we have C =0.267. From (10), (16)— 
(19), and (21), we find 


ri(v) =CMBT~ NGA, T) OK Nat Flo oP K,N,-1} - 


1 Cf, Russell, Mt. Wilson Contr., No. 447; Astrophysical Journal, 75, 340, 1932. 
"2 Handbuch der Astrophysik, 3, Part I, 289, 1930. 
3 The Internal Constitution of the Stars, pp. 323, 335, 1926. 
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When second ionization is negligible, we have 


n(v)=CT 'G(d, T) KN, « (22) 
eee 
From (14) and (15) we find 
_3[{?P 
ugritv)=C.3T * ( Pdp, (23) 
ugt;(v) =C3p= Cyug > N . (24) 
We have also 

P= (p—P) oedorx, : (25) 


In an isothermal atmosphere the K”’s are constant, and the V’s are 
integrals of rational functions of ». When x elements are present, 
these integrals, even when only one ionization has to be considered, 
involve an irreducible function of degree 7 —1 in the denominator. 
It is easy to prove that the roots of this function are all real and lie 
in the intervals between the successive values of K for the various 
elements" so that the integrals lead to logarithmic functions; but the 
general solution is practically useless if ”> 2. 

5. The solution for 7=2 has been given by the writer and by 
Milne.© It may easily be extended to cover the case when a third 
element is present, but not at all ionized. 

If e€ and é’ are the abundances of the two active elements XY and XY’, 
we have 
(P+K)(P+K’) 


p_ TA) TA 
=i+ O(P+R) ’ 


p (26) 


where 
O=eK+e'K’ OR=(e+e')KK’. 
4 Two elements with equal K may be lumped together. 


13 Mt. Wilson Contr., No. 383; Astrophysical Journal, 70, 76, 1929. 


%6 Monthly Notices of the Royal Astronomical Society, 92, 166, 1932. 
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We then find 
ugN =ep, 


\ 


KeP {_ , P+K’| 


rae A Secale P 
ugNi= 0 \1T pyr { (1—e)K log (+7) | is 
27 
ée 6? P\ | 
Total K log (1+) "4 
, P(P he aa : >) ) 
ugNo= o } - ra JLO-k | +(1—e)K log (147) | 
: | (28) 
epee ria (145 
ete’ K 6 a | 


The values of V/, Ni are obtained by interchanging primed and un- 
b] os o oD 
primed letters. 
We now find, after a little reckoning, 


ug(KN,+K'N}) =P?+P(O—K—K’+R)+(1—©)K? log (:+7) 
2s : (29) 
+(1=e)K" log (1+), )—R? log (1+) , ; 
and, by (22), 
m=C.MT*G(A, T)(KNo+K'N}) . (30) 
A direct integration gives 


eigee *s 2P34+3P?(0+K+K'—R) 
0 (31) 


: fs ( P P ] 
+R(K—R)(K —R)} log (+2)-piri| , 


The proportion of inert gas is €¢,— 1—¢e—e’. Let us now suppose that 
a fourth element, X’’, is present in so small a proportion, e’’, that 
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it makes no sensible contribution to the pressure. We then find 
easily 


ug” = ep 
eK" 


yer 
ugn; 0 


| 2P+(0-+K-+K'—R— 2K") log (1+ jen) 


R(K—R)(K’—R) { P(K"—R) P 
~~ (K”—Ry | R@+R) 1° pe | (32) 


+log (+ )}| ; 





When K” =R this goes over into 


i da add : mo / P. ) 
ug; a 2P+(0+K+K’—3R) log (1+) | 


Q | | 
P(P+2R)) | 


+1 BEB) oe as | 


6. When there is but one active element, Y, we may set €’ =o. 
Then O=eK, R=K’, and 


wgN= et (rt)P (34) 
wgN.=2P—(1-6)K log (1+) (35) 
wgNe= pe —(1- 0) P=K log (14%) | (36) 
r1=C:MT #G(A, T)KNo , (37) 
ee ‘ Rt +OP) (38) 
ug =" 2P+ (1-+e)K—2K"} log (x+ z)| | (39) 


” idl tid | P2 ; . ae P P ) , 
ugn, =k Kt (1 +e)K — 2K | glee (+57)} | . (40) 


i 
) 





Equations (34)—(36) have been given by Milne." 
17 [bid., p. 169. 
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7. If X’is completely ionized, we may set K’=~ , R=[(e+e’)/e']K, 
and obtain (after expanding certain logarithms in series) 


P 


, € { a ‘\pD * / ! 
ugn ae” P+R Fe )L +AK+€eR} ’ (41) 
+ ee : ( F , F 
ugNi= PLR (1—e)A log +) )+R log (2+5] : (42) 
e Ff P ) 
T — f / > / ) — , y ny j 
ugn. / PrRiitePreR} +1 e)K log (14%) | s 
43 


—R log (1+) ; 
ug(KNo+K'N,) = (1+) (3P*+5 KP) 
+(1—€)K? log (:47.)-® log (:+7,) ‘ 


tT is then given by (30). We have also 


a 


NU +N =" N, 
€ 


ie P\ , R(K—R) {P(K"—-R) } 
nen ote) log (1+ cn) + (RHR | RFR) 

P rot tie 

—log (x+5,)+log (1+5n)} | +t 


/ 


If there is no inert gas present, e+e’=1. If then the completely 


ionized material is in great excess, our equations become 
ugN =2eP, 
| 
N K | m2 | 
ug, = 2e og (1 ) | 
K (46) 


\) | 
ugN.= a P—K log ( +7 : 


which are those of Milne’s problem III." 


8 Thid., 89, 26, 1928. 
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8. We may now test Milne’s rule that the quantity n= PN,/2KN, 
is not far from unity. He has proved that 7=1 exactly when only 
one element is present. The deviation will be greatest for an element 
present only in traces. In this case we may use equations (39) and 
(40), which assume the presence of two other elements, one partly 
ionized and one inert. 

If we set P/K=y, P/K’=y", these may be written 


ugNy! = : 2P \yfly")+Uutefy’)} , (47) 
yr é” aa { 4; ‘a tI\) 
ugNo = € Kt vhs (y )+Ut+ehly )} ’ (48) 
when 

f(y) =, log (t+y)=3—29+éy ...-, (40) 
ia -S log (r+y)=3—3ytiy...., (50) 
fy) = a oe log (1+-y)=4—dy+3y ..... (51) 

Note that 
sfi(y)=1—-2fily),  swhly)=1—-2f.(y). (52) 


These functions are given in Table II, in which 1/(2+-y) is chosen 
as argument to make the tabular differences run smoothly. When y 
is large—that is, when the active element YX is slightly ionized 
n=f.(y’)/f,(y”) and ranges only from 3 to 1. When y is small, n= 
fily’’)/f.(y’"). This becomes large when y” is great. The error of 
Milne’s approximation is serious only for an element which is just 
beginning to be ionized in the presence of an excess of highly ionized 
atoms. We have by (52) 

_ tet fay—(1t0)y"}fuly”) 


nN T% 


~ i+e+ lay—(r+e)y" Ay”) * 





Milne’s approximation is therefore exact if 2v=(1+.e)y”’. Now by 
(21), y’/y=Bi’ B,/ Bi’ B, e-"?”?, Hence this happens when Y”’ 
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is slightly easier to ionize than XY; but J—J” will be only a fraction 
of a volt. For the active element YX itself, y’’ =y and 


Oi (1—e) f(y) 
1—(1—€)fa(y) ’ 
so that 2 =1 when y is large or small. It has a minimum, 0.93, when 
€=0, ¥=3. 
TABLE II 
VALUES OF f(y) 

















y | 1/(2+y) fr(y) | fa(y) f3(y) 
| 

Oi, ; | ©.00 0.000 (| ©.000 ©.000 
98... 5 OI .023 .O10 .O10 
AB ave | .02 .O4I o1g .020 
£5y | .05 o80 .046 .050 
ae ‘ Fo" 137 .OgI . 102 

Pe ss a | JES .186 tae ey 

ae | 20. 231 .179 214 

> Tan 25 275 | {gas 295 

ie. < . 30 318 | . 273 . 340 

OBO. 5: Be, 301 .324 411 

O50; .. .40 4oo | .378 488 

ee 7 452 | -437 572 

0;003..: | 0.50 0.500 | 0.500 0.667 


Equations (26)-(30) may now be simplified. If we set 
e=a(1—&) , é =a'(1—-&) , (53) 


and P/K=y, P/K’=y', P/R=y,, we find 


/ 


a+a’=1, y,=a'y+ay’', (54) 
—_" 1+y)(1+y’) 
ugnN =al (:-«+! eed , (55) 


ugrs=CMT 'G(d, T)P?{1— (1-0 f(y) — (1 aly +L.) } > (56) 


, y  i1+7 
weNi=P {a(2 427° (37) 


re a’y Kp | 
Vr hy.) 724 Oh) +2 my filyr) , 
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/ _P? { a(1+y’) ee : 
ugNo=_ i+y, (1 e)f2(y) +a’f.(y,) { ) (58) 
* 3 —% D2 ( sail 
peel PY, gy | 
I € Vr | 
(v-—y)(vr—-y’) (1 ) (59) 
< ABE 2 + Laie? 8 — a | 
. Vi ( I +r fy ") ¥ ) f ; } 


g. Consider now an atmosphere in which practically all the opacity 
arises from the interaction of electrons and ions produced from a 
single active element, X, of abundance e—the rest of the gas being 
inert. The number of atoms effective in producing a given line may 
be taken as that of the atoms, in the appropriate state, which lie 
above the level at which 7 has some fixed value 7. Equations (35)- 
(37) may then be written 


ugN = 2P{z—(1—-e)fily)} , (60) 
P? ; 

ugNo= 5 i11—(1—e)f,(y)} , (61) 

ugr.=CMT *G(A, T)P?{1—-(1-fxly)} (62) 


The number M of atoms in any state, with excitation potential E and 
weight g, is given by 


M=5 Ne-BAr , (63) 


Setting 7=7,, eliminating P, and using (21), we have 


M,= Sods L337 e2:- Eo) TGA, 7)" : (64) 
M,=2(ug)* 2 ( a TGA, T)*f(y,€), (65) 
é B, a b] J\7) , \ 2 

where 


fy, )={1-G-O)fi(y) }1-G-O)f) +4 » (66) 


which always lies between 0.707 and 1. N, and N, are obtained by 
setting E=o. For hydrogen we must substitute H(7) for G(A,7). 
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From these equations it follows that: 

a) The abundance of X has no effect on the intensities of its arc 
lines, and only an insignificant direct influence on its enhanced lines. 
The inert atoms, however, may have an indirect effect by altering 
the mean atomic weight uw. Dilution of a metallic atmosphere, too 
cool to ionize hydrogen or helium, with these gases strengthens the 
enhanced lines. The arc lines are unchanged. This effect approaches 
a limit as the dilution becomes great. 

TABLE III 


VALUES OF a(Vz) 


° I 2 3 4 | 5 6 7 
BERS fo ericgvs Bt >|. O.02 0.08 0.22 0.46 0.72 °.90 0.98 I.00 
30. : .02 00 15 35 62 .84 95 °©.99 
40 ae .OI 04 12 29 55 79 g2 0.98 
50 .| 0.01 0.03 0.10 0.24 ©.49 0.74 0.90 0.97 


b) The surface gravity has no effect on the arc lines of X. The en- 
hanced lines show a positive effect (stronger in giants) of normal 
amount (varying as g~') and the same at all temperatures (barring 
higher ionization).'? 

c) The temperature influences several factors. If s=1/vT and 
w=FI,/V, we have 


dlog G(X, T)_ ——_—«V2?— Va-+1—e-¥ ne 
Vds ~ V2? Vse+(Vstwl?22)e- 2 2 . 


When z=0, a=1/2w; when s=«,a=1. For the observable spec- 
trum, V lies between 1.6 and 4. With FJ =8o0 we have to deal with 
values of w ranging from 20 to 50. The corresponding values of a are 
given in Table III. We have then, by (64), 


d 


Ps log M,=vT —aV+I—-E,. (68) 


9 Cf. Milne and Chandrasekhar, op. cit., p. 170. 
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Hence M, decreases as J increases, unless E, +aV >/,+277T. Now 
the excitation potential for the upper state of any line is E+ V, since 
V measures the energy drop in the line itself. This value is less than 
I, except for a few states in which two electrons are excited. All arc 
lines (except the very few having such upper states) will therefore 
steadily increase in intensity with falling temperature, until the 
ionization becomes so small that other sources of opacity have to be 
considered. 

For hydrogen we set d log H(T)/dz= (sz) and find (determining 


d log S/dz graphically) 


s=—p, fe) 0.5 I 2 2 4 6 8 10 
v1 

Boo... O02 O06 0.22 O90 4.83 21.68 3.67 1:63 1:6 

Dang. 3% 0.04 1.61 3.58 3.86 3.71 3.63 3.54 3.50 3.48 


(The values 8’ apply beyond the limit of the Balmer series.) We have 
then 

: = a ee a (69) 
For cases in which most of the electrons come from hydrogen, the 
lines of the Balmer series should therefore decrease steadily with ris- 
ing temperature. 

Milne finds a minimum for these lines because his formula for K 
contains 7~! instead of 7~? (at lower temperatures, T~ “'). This 
substitutes 7? (or T3) for T~ in the expression for M,, which then 
still increases at high temperatures. The weakening of the hydrogen 
lines in classes B and O can be explained on this theory only if there 
is a large increase of opacity due to the ionization of helium and 
higher ionizations of other elements, a condition which demands 
that these elements shall be more abundant than the hydrogen. 

For the enhanced lines, we have from (60) and (49) 


{ I ) 
d log N,_ d log e* Ay Aa 2+ay | d log y 


dz dz 1—(1—e)fi(y) dz ’ 
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to 
wn 
U1 


and from (62) and (50) 


\ 


I 
d log P ‘i sities ae a ) dlogy_ 


8a ie 3, — 
ieee 1— (1-6) fly) ds 7° 


also 


>) 4 ay 
dlogy_dlog P_dlog K_d log Pa ra sor. 
dz dz dz dz ° 


We then find 


(1-—e) /fily)—f.(y)} 


d log N, | 3 
~~ os 1—(1—e)fi(y) 


7 30T —LaV+ (I+3eT—JaV), 


and may write 


d log M . a . gist ; 

a ‘= —E,—30T—3aV+y(y, )\(U+50T—falV) . (70) 
The factor y is always small. It diminishes the rate of increase of the 
enhanced lines with rising temperature, but will not reverse it for 
any known lines (see §15). 


TABLE IV 
VALUES OF 7¥(y,¢) 





| 1/(2+y) 
€ | ae J ae ee 
| 0.00 °o.10 o.20 0.30 0.40 0.50 
fe) | 0.000 | 0.053 | 0.068 | 0.066 | 0.047 . 000 
re 0.000 | 0.025 | 0.030 0.027 0.017 000 
| 
Second ionization, of course, causes all enhanced lines to have 
maxima. 


1o. Consider next a constituent Y’’ which is present in such 
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amount that it contributes practically nothing to the electron pres- 
sure. From (47), (48), (60), (61), we have 


Ny _e’" K hv") +tO+tefi(y”) ' 
N, e K” 1—(1—e)f,(y) F is 
Nye” yf.(v")+ (1 +6) f(y’) (72) 

. le 


Ni € 1—(1—e)fi(y) 


The intensity of the lines of X’”’ depends on its abundance, €’’/e, 
b / 
relative to the opacity-producing element Y. Call this a’’. We may 
now distinguish four typical cases. 
Case 1.—X and X” both slightly ionized. Using (52) and Ky= 
fon) ‘a oD . ¢ 
K’’y”, we find 
N,a"'K" 


Nv=Noa", 9 NY= = Nae (73) 


The arc lines are unaffected by gravity; the enhanced lines show a 
normal positive influence. For the temperature effect we find 


d 


Zc log Mi’ =vT+I-—aV—-E,’ , (74) 
d ” Set 3 ? ” Ad 

F log M;’=—f0T —gaV+1—-I'—-E; , (75) 
az 


the term in y vanishing, for both small and large y. The arc lines 
strengthen with falling temperature, unless /j’>/+v7T—alV. 

If 7” is considerably less than J, the ultimate enhanced lines grow 
stronger as the temperature falls, and so may subordinate lines for 
which £7’ is small. This increase has a limit. When the temperature 
falls so low that NV{’ becomes comparable with V,, X”’ takes over the 
role of active element. At lower temperatures its enhanced lines 
weaken. 

Case 2.——X slightly, Y” heavily ionized. We have 


rer T ‘? 
ivi =3N,a y. 


kK" y; 


NY =3N,a"" 
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But by (60), (61), y= P/K =2N./N;. Hence 
7 K No rr WIN / 
N = 4a” K"” N, ) N; =a'N O- (76) 


The enhanced lines of X’’ now behave in all respects like the arc lines 
of X. They are insensitive to gravity and weaken rapidly with rising 
temperature on account of the increasing opacity due to the ioniza- 
tion of XY. For the arc lines V, « (ug)?. These lines show a negative 
absolute-magnitude effect. Their decrease with rising temperature is 
unusually rapid. 

Case 3.—X heavily, X” slightly ionized. 


2 2 NE par: 
NO =" gory HUN (77) 
yr at lo ( 4 Pa Kk" Ni log ( ie /") 
ae ENE LS ee ee 
Le Gs 
=e log (r+y”) , (78) 


by (60), (61). Here the arc lines behave as enhanced lines usually do 
with regard to both temperature and gravity. (This happens for the 
Balmer series in the cooler stars, where only the metals are ionized.) 
For the enhanced lines the absolute-magnitude effect is twice as 
great as usual, and the increase with temperature is rapid. The op- 
tical depth 7, has disappeared from (78), although it enters indirectly, 
since y’’«7,'. Under the assumed conditions, X’’ is considerably 
ionized only at high levels, and a considerable increase in 7, takes in 
many more neutral atoms, but few additional ions. 
Case 4.—X and X” both heavily ionized. 


N.a"'K 


Kk"? Ni’ =N,a" : (79) 


Ath — 

Ni" = 

The arc lines are again insensitive to gravity and the enhanced 

lines normally sensitive. The temperature effects are the same as if 
X”’ were the active element. 
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Equations (71) and (72) show that the transitions between one 
and another of these cases will be smooth. They suffice for exact cal- 
culation in any desired case. 

Equations (73), (76), (77), (78), and (79) depend only on the as- 
sumption of an isothermal atmosphere and are true for all laws of 
opacity. The dependence on gravity and temperature will differ 
with the assumed law; but the reversed sensitiveness to gravity for 
arc lines in case 2 and the increased effect for enhanced lines in case 
3 are general properties as long as the opacity depends on ionization. 
These effects have been previously described by Milne,” who has 
shown that the relations are qualitatively true, though to a less de- 
gree, when X”’ and X are of comparable abundance. 

t1. Second and higher ionizations remain to be considered. Under 
stellar conditions the second ionization of any element is insensible 
until the first is practically complete, and, when it happens, there is 
always a large excess of electrons arising from other sources. Equa- 
tions (77)-(79) may therefore be applied in all cases, substituting 
I’ for I’, NY’ for Nv’, and NY’ for Ni’. We may distinguish three 
groups of elements: 

i) I’’<13.6. Cat, Sr*, Bat, and, to a less degree, Sct, Y*, 
La*.—Second ionization precedes that of hydrogen but begins 
when the first ionization of the metals is practically complete. At 
this stage we have case 3. Ultimate lines of X’’+ are normally sensi- 
tive to gravity and almost insensitive to temperature. Those of 
X’’++ (if observable) would show exceptional sensitiveness to grav- 
ity. If the ionization of X* is far advanced before that of hydrogen 
begins, we have case 4. The lines of X’’+ are now insensitive to grav- 
ity, and decrease with rising temperature. When hydrogen begins to 
contribute more electrons than the metals, we pass to case 2. The 
lines of X’’+ are stronger in dwarfs than in giants and decrease very 
rapidly with rising temperature. Finally, when most of the hydrogen 
is ionized, we return to case 4. 

ii) /’’=13.6 (approximately). 7i+, Zr+, and the rare-earths. 
Lines of X”’+ behave as hydrogen lines of the same excitation poten- 
tial would do. Starting with case 3, before the hydrogen is ionized, 

20 Monthly Notices of the Royal Astronomical Society, 89, 166, 1928 (case 3); ibid., 92, 
161, 1931 (case 2). 
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they fade rapidly as ionization goes on and the opacity increases, and 
become insensitive to gravity. 

iii) J’ >13.6. Mgt, Sct, Fe*+, the light-elements, and all higher 
ionizations.—-Second ionization begins with case 3 and goes over into 
case 4. Ultimate lines of Y’’+ at their maximum intensity are nor- 
mally sensitive to gravity, but lose this as they fade out; subordinate 
lines lose most of the gravity effect before they reach their maxima. 
Ultimate lines of X’’++ would be exceptionally sensitive to gravity 
at their first appearance. By the time that the highly excited sub- 
ordinate lines (which alone are observable) appear, this sensitiveness 
drops to about normal. When lines of X’’* and X’’++ are both visi- 
ble, they will all be stronger in giants than in dwarfs, but the effect 
will be more pronounced for the higher ionization. 

12. The observable lines of neutral helium will behave in the same 
way, being sensitive to gravity at temperatures below their maxima, 
and insensitive as they fade out in the hottest stars (cf. a Cygni!). 

For the Balmer lines we have / = 13.5, E=10.1. If the atmosphere 
were composed entirely of hydrogen, they would increase steadily, 
down to the lowest temperatures; (69) gives 


d 


1 log M=3.4+vT7 —B=-+5 (for B stars). 
dz 


With metals present, this increase will cease as soon as they contrib- 
ute about half the electrons. When s has increased by o.2 (or T 
dropped by some 1500°), the hydrogen will be almost all neutral, we 
shall have case 3, and by (70) and (77), neglecting y, 


d 


/ log M = — 301 —}alV —10.1= —10 (classes F and G). 
dz 


When the temperature has fallen so low that most of the metallic 
atoms are neutral, we have case 1 and, by (74), 


: log M=0T—aV—10.1+1=—6, 
des 


setting J =6. 
The intensity of the Balmer lines will therefore have a rather 
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sharp maximum near the point where ionized hydrogen and ionized 
metallic atoms are equally numerous, from which it will diminish 
slowly with rising temperature, and much more rapidly with falling. 
| At low temperatures the rate of diminution (for equal changes in 
| 1/T, or the color index) will fall off. The lines will be sensitive to 
gravity in the region where the intensity falls most steeply, but not 
elsewhere. 

The presence of several metallic constituents will smear out the 
transition between the second and third portions of the curve with- 
out altering its general shape. 

13. The existence of a maximum intensity for any stellar line de- 
pends essentially on the presence of other sources of opacity than the 
ionization of the atoms or ions which produce it. This circumstance 
is implied by Milne’s equations, but deserves emphasis on account of 
its importance. 

The general problem of a maximum, where two or more elements 
are being ionized, is complex; but the case in which one element_X is 
partly ionized, in the presence of any amount of material which is ) 
completely or insensibly ionized, is simple. In this case y’ =o in (54) | 
and (59), and we have 





( 
iVN= “yd ad —_ \_ “ay 
ugN =aky hie eo akyu, (80) | 
ugNi=Ky | mre 2f,(a’y) — 2(a’ + ae) f:(y) =Kyu,, (81) | 
| 
oe a ee ee oe = 
ugN,=Ky pep OOO y)—(a Faeo)faly) | =Kyu,, (82) 
ia = K*y?{3(1+a’ —a’e,) +f.(a’y) — (a’ + ae.) f.(y)} =K2y’u., (83) 
CMT *G(A, T) 
Mei Key | Hata! a'e) +uf(a'y)= 27} Kya. (84) 
cw \ i+ey 
Note that 


U2=Ut+U; . (85) 
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We may take y and ¢ as independent variables. By (49) and (50) 


Ua MN = spam i Ala’ =, 
whence 
d(yu)=(1+y)Y dy, d(yu;)=a Vdy, 
d(y?t) =ay Y dy, d(y*u2)=y(1+a'y) Ydy, 
d(y*u3)=a'y(1+y) Vdy, 
where 


(r+ y)(1+a’y)?V =1+2y+a'y?+(1-6)(1+a’y)?. 


We now have 


d log M,= —(I+ E+ 307) dz ae “) —d log g, (86) 


d log t1= — (I+ eT av as) — dlog g=o, (87) 


whence 


dlog M,= {U.(2I+ jxoT —aV) —(I+E+$2T)\dz—(1—U;)dlogg, (88) 


where 
’ au, 


U.(y, a’, € . 
2() ) = (1+a’y) ut 


The numerator and denominator are linear functions of €, and the 
latter is finite; hence the values of U, for « between o and 1 are inter- 
mediate between its values at these limits. 

If a’ becomes very small, while a’y remains finite, y is great and 
f.(y) negligible, so that w= 1/(1+a'y),u.=3+/f,2(a’y), and U,(y,0, €) 
=4+f,(a’y). Ifa’ is small and y finite, V,=1. When a is small, 


we have 


7h, laa 


I (259? rs 
Iyay 


ra / | u 
fly) —a°f.(a'y’) =- yy? log 1+a’y 


y 


Uo = a(2—&)f2(¥) , Uy = 2a(2—e) f(y) , Uz=1I— 36, (89) 


I 


U.(y, 1, €) = 2(1+y)f2(y) ’ 
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Values of U, are given in Table V. These quantities change remark- 
ably little with e, and a’; they depend almost entirely on a’y, that is, 
on the relative number of electrons liberated by the active and the 


TABLE V 


FOR MAXIMA OF SUBORDINATE LINES 


1/(2+a’y) | U2(y, €, 0) | Ualy, 4, 0) | U2(y, 3, 1) | U2 (y, €0, 1) 
0.00. 0.500 | 0.500 =| 0.500 | 0.500 
02 190 | 526 | 530 | 533 
05. 540 | an 563 | . 566 
10. 591 | 602 604 | .613 
BS. oss: | 645 | 647 | .656 
20.. 679 .690 695 .698 
25. 725 | Sa 739 | .740 
30... M4 | ‘779 | 785 | - 783 
as. .82 | 27 | .830 | 830 
AO. . .878 | 881 | .884 882 
BSc 0.930 | 0.9360 | 0.937 0.937 
0.50.. 1.000 | 1.000 | 1.000 | I .000 


wholly ionized element. This ratio is N,/N;=w,/a'u. When a’ is 
small, 


N 2 a’y i(a’y 
: oe (rta'y)fi(a’y) ba 


I 
p a’y 
As this ratio increases, the sensitiveness to gravity diminishes, some 


corresponding values being 


N, 

Ninn Q.00 0:22 0.43 1.05 1.56 4.4 oO 
4¥1 

1—U,... 0.50 0.45 0.4! 0.32 0.23 0.12 0.00 


When a’y is large, d log M/dz is negative; when a’y is small, 
d log M/dz is positive, unless E+aV >/+v7. This condition is 
satisfied only when the atom in the upper state contains more 
than enough energy to ionize it. All other lines will have maxima 
when 

[+E4 301 


2I+40T —aV ' (91) 


U,(9, a’, €) = 
If E is small, these arc lines will still show a pronounced absolute- 
magnitude efiect at and beyond their maxima; if F is large, the 
etfect will be small. 
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14. For the Balmer series, we may take V = 2.8 (near Hy), T= 
10,000, whence 77’ =0.862, aV =1.18, U,=0.894. From Table V, 
setting €,=a’=o (since both are known to be small), we find 
1/(2+a’'y) =0.414, a’y=0.42, and, by (go), 


N,_ 2.18 
Ni a’y 


$43 


the hydrogen, therefore, produces five times as many electrons as the 
metals. 

The assumed temperature is rather high. If T= 7000°, U,=0.932, 
a'y=0.245, N,/Ni=8.6. 

Since hydrogen furnishes most of the ions, it is better to substitute 
B for aV in (91). This gives T=10,000°, 7000°; a’y=0.56, 0.48; 
V,/Ni=4.0, 4.6. There is strong evidence that in stellar atmos- 
pheres a’ is 0.01 or less. The maximum of the Balmer lines, therefore, 
comes when the hydrogen is not more than 5 per cent ionized. 

For any other element we have 


yl’ B, f qo 
y ~~ 2B,e!"—Ds (since 2 for 1) 


and, approximately, Vj’/N(’=2/y’. For iron, /’’=7.83 and B,/B, 
=0.75 (at 7000°). With the Balmer lines at maximum at 7000", 
Nv/Nt/ =15,000 a’, while at 10,000° it is 1000 a’. Hence, unless a’ 
is extremely small, iron and all other common metals will be very 
highly ionized, and our basic assumption is verified. 

For the higher ionizations in B stars, hydrogen is in great excess, 
and we may take ¢,=o, a’=1, and neglect 8, which is small. We 
now have V,/N.= 2f,(y)/yvf.(y), and find by computation the accom- 


I I 1 U2 y N;/No 
He 4471. | | 2h 20.9 | 1.5 0.906 | 0.37 | 5.7 
Ile* 4686 Sh ot. Gieea 48.2 2.5 0.927 0.27 a5 





panying tabulation. For maxima of observable lines the ionization 
is always high. This is Milne’s problem III (modified by the change 


in the opacity formula). 
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15. For enhanced lines of Y we have 


d log M,=—(1+E+ gor yds 20") —d log g. 
The condition for a maximum is 
=... — eA _y , E+ioT+jaV : 
(t+a’y)ur Us(y, a, 0) = 3+ 2I+40T—aV ° (92) 


When a is small, we find, from (89), 1/U,=4(1+y)fi(y); U, ranges 
from 3 to o as y goes from 0 to «, and there is no real solution. 
When a’ is small, we have 
_ FA fla'y) ef.) 

1+ 2(1+a’y) }fi(a'y) —eofs(y) } ° 





U; 


Now f(y) >f.(y). If a’y is finite, f:(y) and f.(y) are negligible and 
U,<%, there is no solution. But if y is finite, a’y=o and U,= 
[1 —eof2(v)|/[2—2€f:(y)]. The maximum value is 0.534 when €,=o0, 
y=2. There will be two real solutions, giving a maximum and a 
minimum close together, if £,+0.63v7+0.53aV <0.034/; but this 
condition is not met by any known enhanced line. 

Since 1—U, is usually greater than 0.5, these lines are unusually 
sensitive to gravity. 

16. The lines of other elements may have maxima and minima, 
owing to the change in opacity, produced by the ionization of X, 
even though they are ionized completely or not at all. In the first 
case N,’, and in the second case N;’, is proportional to NV, the 
“amount of matter above the photosphere.” If M” is the number of 
atoms of either sort which are excited with potential /’’, we have 
from (80) 

d log M" =—(I+E"+ gor yds 90" a log g. (93) 
Eliminating y by (87), we find 
d log M={U1'(2I + 40T —aV) —(I+E"+$0T) }dz—(1-—U"’)dg, (94) 


when 
Hy at @\ ae AEE) ta 
a at cere 
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Here again U2’ varies monotonically with €,, and 


UY’ (y, a’, 1) =3+f.(a’y)-—f.(y) , (95) 


_Gty) 1+a’+ 2f,(a’y) — 2a’f.(y) } 


Ut’ (y, a’, 0 
2(y, @, 0) 4+2y+2a’y 


(96) 
The first of these quantities is given in Table VI, and also the 
difference of the two. This difference is small, and the presence of 
inert gas again produces very little effect. 
For a stationary value we have 


yin THE + or (07) 
oat — oF 
If this has any solution, it has two. The larger value of y corresponds 
to a maximum of M”’, the smaller to a minimum. Between the two, 
the increasing opacity due to the ionization of X “raises the photo- 
sphere” and diminishes M’’. If a’ is large, there will be no such 
effect, even for ultimate lines; but if a’ is small, it will appear even 
for considerable values of E”’, while for ultimate lines the maximum 
and minimum will be widely separated. 
When a’ is small, we have, approximately, 


Ul = (r+y){3+f2(a’y) —ef2(y)! 
s y+(2—6)(1+a’y) 


If y is finite, this becomes 


(y+1) {1 —Gf2(y)} 


Uy = 
VF 2a--& ‘ 


(98) 


but if a’y is finite, we have 
U! = s+f.(a’y) : (99) 
Each equation gives an approximate value of one of the roots. In 


the interval in which y is large and a’y small, U;’ is very nearly 
unity; hence practically all the lines show maxima and minima. 
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17. If the opacity arises from interaction of electrons and neutral] 
atoms, we have 


se 


d log 72= — (21+ 307) dz pe: —d log g=o, (100) 


d log M,=|U;(21+ 30T) — + E+ joT) |dz—(1—U;)d log g, (101) 


where 
au; (1+a’y)U,—a 

Ta(ity)u a’ i+y) 8) 
When y=o, U,=1,and when y= ~, U,=4,so that there are maxima 
for practically all arc lines, as before. When a’ is very small, U, 
ranges from 3 to 3 as a’y decreases from large to small values, and 
then from 3 to 1 as y does the same. If a’=o, the first part of this 
range is lost, and there is a maximum only if E>3/—§v7. Other- 
wise, M increases steadily with increasing zs or falling temperature. 

For an element of unchanged ionization we have, by (93) and 
(87), 


log M" = 1U3/(21+ 307) — +E" +307) | dz—(1—U7/)d log g, (103) 


where 
rer 3 
U3 ~ alu" (104) 
When y is small, 
eo = ay 
™ " 12— 0¢€ 


and when y is large, 


l rer = 1 
2a'(1 +a’ —de)v- 
As for the other law of opacity, there will be maxima and minima, 
provided that E does not exceed a certain limit, increasing with a’. 
When a’ is small, we have 
on _tta’y}{3+f(a’y)i—1_,, 
U;= =U,. (105) 


/ 
ay 
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As a’y decreases from large to small values, U{’ goes from 3} to 2; as 
y decreases to small values, U{’ returns to 3. Hence, if E<}/J— v7, 
there is at higher temperature a maximum in the first range and a 
minimum in the second; but if & exceeds this limit, M increases 
with 7 (or as z decreases) until opacity due to the ions becomes im- 
portant. 

When both types of opacity are sensible, as may occur in a slightly 
ionized hydrogen atmosphere, we have 


ug(ti +72) =C ST : K*y?(Hu,+ Aus) , (106) 


where A =C,/C,a’(1— 6). 
Since dH /dz= 6, we have 


d log (t:1-+72) = —(2I1+40T—<x8) dz | 


H(1+a'y)+Aa'(i1+y) Ydy » (107) 


T Hu.+ Au, HAO Ss 


where x = Hu,/(Hu,+Au;). Then 


dlog M,=U.,;(21 + $0T —xB) —-(U+E+ $07 )dz—(1—U,;)dlogg, (108) 


if 
U.- H(ita’y)U,+ Aa’ (ity)U, 
8 W(r+a’y)+Aa’(it+y) 
If 
& — H(1+a'y) moe 
1—& Aa'(it+y)’ ae 
—** : iU, 
U.3,=§U.+(1—§)U;, a : (110) 
For X” we may use (108), substituting for U%,, 
Us =tUY+—pur= UTM y,.,. (111) 


au 
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From (12) and (14), Aa’/H=k,/(1—€.)k,, Where x, and x, measure 
the opacities due to equal numbers of ions and neutral atoms. Set 


3 2 G; 
inne, 112 
gy Caney — 
This ratio varies with the temperature but is probably always small. 
We then have 


" ee 
as (1 —€)) U2; Us; } 
a } 


Uy= uy b+1 —e+(b+a'—a’e)y\ , (113) 


For a hydrogen atmosphere a’ is very small. When y is finite, we 
have : 


Uo = U2 = 1 — €of2(y) , u;=a'(1—e+3y) , u=y+2—-6. 


Taking €,=0, we find 
r+b+2by 


1+b+by (114) 


U,,=U;= 


= 1+b+%by. 


For the maximum, 


i (1+6)(/—E+27T— xB) ) 
by =~. — , 115 
“ E~3I+G0T +548 : 
When £>3/ —{v7'— 3x8, the maximum occurs for a finite value of 
by, independent of a’ (to the first approximation). But if we suppose 
that a’y is finite, we have 
I 1 r , / 
u,=otf2(a y) , U;=U2— Uo , u=y(1+a’y) , 


Ita’y ; 7 


Uo = 


and find, when €,=0, 


(b+a’)(1+a'y) |3+f.(a'y)| —b 


a’ '\1+b+(b+a’)y!} (116) 


U,=Us= 


and 


I b I 


ma g--— 3 


(I +a'y) 141,(a'y) . 
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When a’/b is small, (116) reduces to (105) and x is negligible. The 
maxima for small values of E come when a’y is finite and are inde- 
pendent of b. Hence, if a’ << <6, the maxima of different lines of the 
same element may be widely separated. When E=3/, the maxi- 
mum will be very flat. 

18. In applying these formulae to actual stars, it is customary to 
assume a constant value of g; but it is very little trouble to allow for 
the variation of g along the main sequence and the giant series. 

For the former, the best values of g are obtainable from eclipsing 
variables. All spectral types from B to M are available.’ For the 


TABLE VII 





| | | | 
| log | A : . | log | : 
Star | Sp. | 2 | (e/t0) | i= 4 | Star | Sp. 2 | ue O— 
} | | || i | 
V Pup.. Br | 0,53 |—0.54|—0 o4|| 8 Aur..... | Ao I.16]/—0.52|—0. 20 
u Her.. B3 | .64|— .42\+ -0s]| TX Her....| Aa 1.3 |— .o4/+ .24 
o Aql B3 | -64|— .gorr .17] 2Her......| Fg 1.7 |— .29\— .11 
Z V¥ui.. B3 64 |— .51/— -o4]| Sun........ | dGo | 2.0 .00|}+ .09 
U Oph B8 go j— .98i+ .12] W OUMa....1. dGo| 2.3 |+ .16— .20 
RS Vul. Bo .go |— .60/— .20]] aGemC...| dMo| 3.3 |+0.19/—0.08 
TV Cas Bg | 0.97 |—0.52)/—0.14 | | | 


brighter components of the best-determined stars the data are in 
Table VII. The formula 
§ 
logo | = —0.65+0.282 (1174) 
represents the data with an average residual of only +0.12 for the 
individual stars. 

For the giant stars we may use Strémberg’s” recent determina- 
tions of the mean absolute magnitude. From his Tables II and II] 
we find the data in Table VIII, for which the average residual given 
by 


=+0.76—0.972 (1176) 


0 


logro 
ul 


re 


g 

g 

is to.07. For the natural logarithms, 
D 


log g= —1.5+0.62 (dwarfs) = 1.7 — 2.23 (giants). (118) 


2" Russell, Dugan, and Stewart, Astronomy, 2, 719, 1927. 


22 Mt. Wilson Contr., No. 442; Astrophysical Journal, 75, 115, 1932. 
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These formulae, although purely empirical, should give a trust- 
worthy representation of the variation in gravity along the two prin- 
cipal sequences of stars. Strémberg’s data for the giants represent 
the stars brighter than the sixth magnitude, which is just the selec- 
tion we need. 

If now we set d log g= Bdz, equation (88) becomes 


d log Mo 


ee =U,(2I+ 3 oT —aV+B)—(I+E+3 oT+B). 


If E+j0T + 3aV <3B, U,<} and there is no maximum——M, in- 
creasing to the lowest temperatures at which Y’ remains completely 


TABLE VIII 





| 

| 

| (D/Do) (m/ mo) (g/g) z cial 
= | ——— 
Ao... -+| 11,200 | +0.6: | 0.27 0.48 |—0.05 1.03 |+0.18 
Ag gooo | +o0.8 | 0.48 45 |—0.51 I.29 |— .02 
gk 3 zi 6000 | +1.5 | 0.56 36 |—0.76 1.78 |+ .20 
gG3 ee 5000 | +0.3 1.03 ga: [553 2.32 |+ .04 
gK1 a 4o0o | —0.2 | 1.32 | 60 |—2.04 | 2.90 |+ .o1 
gks --| 3300 | —1.3 | 1.71 | 77, |—2.65 3.52 |? .02 
gM1 ae 3 3000 | —2.0 1.93 | 0.92 |—2.94 3.86 |+0.05 

| | | 


ionized. This state of affairs does not happen when the absorption 
arises from metallic ions, but it does happen when the absorption 
comes from hydrogen ions or neutral atoms, and aV is replaced by 
B or by o. 

19. We may now apply our formulae to the maxima of lines in a 
hydrogen atmosphere containing small amounts of ionized metals. 
Here €,=o and a’ is small, so that we have U,=U!/ =}+f,(a’y) and 


jy E+oT +484 4B 
flay) =r ter B+B (119) 


When the metals produce most of the ions, we must substitute aV 
for B. 

Assuming as before V = 2.8 (A 4400), we obtain Table IX, which 
gives data for giants and main-sequence stars and for the assumed 
temperatures at which the maximum occurs. The upper part of the 
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table gives the values of a’y derived from (119), which indicate the 
degree of ionization of the hydrogen, and of a’y/’ and a’y;’, which de- 
fine the extent of the first and second ionizations of the metals. For 
the maxima of 77+ and Fet, a’y>2; most of the ions come from the 
metals, and aV has been used in (119). For JJ and Mg* the reverse 

is the case. 

TABLE IX 
MAXIMA AND MINIMA FOR LINES IN HYDROGEN ATMOSPHERE; 
OPACITY DUE TO IONS ONLY 


Element H Tit | Mgt Fe 
I;. | 13.54 | 6.80 | 7.61 7.83 
I; tee 13.6 14.07 16.5 
E 10.15 r.26 8.83 2.80 


gee 10,000} 7000 10,000 7000 10,000 | 7000 10,000 7000 


Maxima, giants: 


a’y 0.63 | 0.55 17.9 12.5 0.98 °.89 6 3 5.0 

a’y,’ ; wee] 2.2XIO—3] 4.9K 10-5} 3.8X 10-4] 1.8X10 2.0XI0 9.6X10-s 

a’ys’ 3 17 5.1 9-5 193 770 

_ } 

Main sequence: | } 

a’y.. | 0.55 | 0.47 8.0 | 6.9 0.84 0.76 4.1 | 3.4 

a’y;’ es 9.6X10-4] 2.7X10-5| 3.3 KX 10-4} 1.5 X10 I.3X10-3] 6.5 X10-S 

a’ys’. rey ene Are 9 1.4 7.6 126 520 
Minima, giants: } } 

y.. 3 , 0.13 «| 3.2 O.44 

eS lass Mens 0.17 16.8 | 13.7 
Main sequence: | 

a, 0.58 : 3.8 | 0.69 

a . 0.70 19.5 2f..4 


We have, approximately, Nj’/ Ny’ =a’y’’/2a’. Hence, unless a’ is 
very small, the proportion of neutral atoms is small, and our method 
is valid. Nevertheless, if a’=0.o1 (about the largest permissible 
value), 1 part in 200 of the iron is neutral, if the enhanced lines reach 
their maxima at 7000’, and 1 part in 10, if the maximum comes at 
10,000°. The arc lines of iron are actually nearly as strong as the en- | 
hanced lines in classes F'5 and Fo, when the latter reach their maxi- 
mum and still appear in large numbers in class A. As the observable 
enhanced lines correspond to transitions of small probability, this | 
relation of intensities is what might be expected. Second ioniza- 
tion, even at the higher temperature, is negligible. 

By (98), minima of these enhanced lines would occur, if second 
ionization did not intervene, when 
2E"+30T+B+B 

I—E”+0T-B 


y= 
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The corresponding values of y and y” are given in the lower part of 
Table IX. Since the second ionization of 77 is far advanced before 
this condition is reached, the lines of 77z* will actually fade rapidly 
as the temperature rises. Since for Mg and Fe it is sensible, but 
small, there should be an inflection in the curve, if not an actual 
minimum. 

We may now compute the difference of the temperatures at which 
the minima for the various lines occur. From (87), when a’ is small, 
we have 

2+a'y 


Gienitia: d log y—(21+ 307 —aV)dz—d log g=o. 


The coefficient of d log y is always between 2.00 and 1.84. For hy- 
drogen, with the assumed conditions, we have, approximately, 


2.9 As=1.9 A log y—A log g=1.9 A log y—B Az. 


If the maximum for Ti+ comes at 7ooo° for giants, a’y=12.5. For 
the maximum of hydrogen, a’y =0.6, almost independent of tempera- 


ture. Hence Aln y= —3.0, As=—o.21. ‘A second approximation 
gives a’y=o0.58, Az=—o.22. Similarly, for Mgt, As=—o0.03; for 
Fe*, As= —o0.16. For the main sequence the differences are 77", 


—o.17; Mg*, —0.02; Fe+, —o.12. For giant stars s increases by 
0.55 per spectral class (from B to K). Hence the computed values 
of Az indicate that the maxima of /e*+ and 77* should occur 0.3 and 
0.4 classes later than those of H and Mg*. The observed maxima 
are Ti+, F5—Go; Fet, F5; Mgt, A2; H, Ao—in the same order, but 
much farther apart. 

Part of this discrepancy may arise because almost half of Miss 
Payne’s® F5 and G stars are super-giants. If gis one-tenth as great 
for these as for ordinary giants, Az would be changed by —o.08. 
Again, the spectral energy-curve of A stars is disturbed by hydrogen 
absorption being relatively low in the red, beyond the limit of the 
Paschen series, and high in the violet—-so that the color index gives 
too high a temperature. McCrea’s calculations™ show that, for stars 
with atmospheres of pure hydrogen and temperatures of 8000° and 


3 Stellar Atmospheres, pp. 127-130, 1925. 


4 Op. cit., p. 855. 
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10,000°, the slope of the energy-curve in the visible spectrum cor- 
responds approximately to temperatures of 12,000° and 10,000°, so 
that this may account for the whole discrepancy. 

20. If, however, the neutral hydrogen atoms produce much more 
opacity than the metallic ions, we must use (111) for H and Mg*, and 
(105) for Fe+ and Tit, with the results given in Table X. We have 


TABLE X 


MAXIMA IN HYDROGEN ATMOSPHERE; OPACITY DUE 
TO NEUTRAL ATOMS AND IONS 


| 
| 











Element | H | Met | Ti* | Fet 
eae —— or — | ; | : | | 
| Ree 10,000} 7000 10,000] SOOO Hl. ia tie o8 10,000] 7000 10,000) 7000 
Giants: | | | | Giants: — | | 
U3.......] 0.855] 0.865] 0.806) 0.814 bt CRN kee. 525] 0.519) 0.585) 0. 580 
by 0.77 |} 0.73 | 1.40 | 2.26 a'y 14.4 |20.2 2.20 | 2.60 
WAN ie it Bee 2.33 | 2.48 | 2.46 hae -....| 2.03 | 2.02 | 2.26 | 2.20 
Ee, i | . | | 
cS } ‘ 
Main se- | Main se- | | } 
quence: | | quence: | | | 
Wigs aie ein | 0.868] 0.879) 0.823] 0.833 | ae 0.569) 0.564 0 622) 0.619 
by/(1+b)..| 0.66 | 0.57 | 1.14 | 1.00 Dien 3.40 | 3.88 | 0.76 | 0.82 
| 1 
now, from (107), setting €=o, 
d log (71 +7.) = — (27+ $0T —xB)dz—d log g +Zd log y, 
where 
gat Viitbt(bta’)y| | 


a’u,+bu; 
If a’ is small, and y and 0 finite, 


7 ltt btby) , | 
“~~ 7+b+2 by ’ 


but if a’y is finite, and a’/b’ small, 


2 +a’y 


ak (r+a’y)Ut’’ | 
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where U%’ is given by (104). Hence Z=2 when by is small, rises to 3 
when by is large and a’y small, and drops to 2 when a’y is small. Its 
values at the maxima are given in Table X. 

Comparing the maxima of H and Mg*, we have (closely enough) 
for giants — 27 Azs+2.4 Alog y=o, whence Az= —o0.05. That is, the 
maximum for H comes o.1 spectral class earlier than for Mgt. 


Similarly for Fe+ and Ti+ (near 7o0o0°), Az = —o0.16, the maximum 
for Fe+ coming 0.3 class earlier., Comparing Fet+ (at 7o00°) and H, 
we have Aln y= —1.1+In (a’/b). In this case Z rises to 3 in the 
interval, so that we may adopt a mean value of 2.7 and obtain 
Azs= —o.11+0.1 In (a’/b). The observed separation is about 13 spec- 
tral classes; Az = —0.8, which gives a’ =0.001b. With a’ =0.01b, we 
would have Az= —o0.57. In view of the corrections already men- 


tioned and the uncertainty of the observed maxima, no definite de- 
termination of a’/b can be made. 

21. These results may be compared with those predicted for an 
atmosphere composed of some element entirely in the first state of 
ionization (e.g., Na) with traces of H, Fe*, etc. In this case we have, 
from (46), 

ugn,=2eK {y”’—log (1+ y")} , 


and from (62), with e=1, 
ugr=CiS3T—? FIK*y? , 


in which we have set G =o, since there are supposed to be no neu- 
tral atoms of the active element. 
For the maximum we find easily 


EE” +40T 


+E" + por POLO"). 


This is equivalent to the equation (38) used by Milne and Chan- 
drasekhar’> in discussing maxima. 

Only the atomic constants for the element under consideration 
appear, for obvious reasons. The lines have maxima, but no minima. 
In Table XI, the entries headed y show the ionization of hydrogen 
when the other lines are at their maxima. The values of Az (com- 


25 Op. cit., pp. 180-183. 
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puted as before) show that, in a once-ionized atmosphere, the maxi- 
ma for Fe+ should come at the same temperature as for Mg*, and 
higher than for H. This is contrary to fact. A still more serious dis- 
cordance is indicated by the values of y;’. At the observed tempera- 
ture of maximum for Fe*+, only one iron atom in a million should be 
neutral. As a matter of fact, the arc lines of the metals remain as 
high as class A. In the spectrum of Sirius, 248 arc lines of various 
metals have been observed, and only 76 enhanced lines.”° This is con- 
clusive evidence of the great abundance of hydrogen. 

The computed maxima for Fe+ and 77* occur when second ioniza- 
tion is moderate, and the maximum for Mgt when it is large, on ac- 
count of the high value of EF. 


TABLE XI 


MAXIMA OF LINES IN SODIUM ATMOSPHERE 


i Ti" Msg" Fe* 
T....|10,000) 7000 10,000 7000 | 10,000 7000 | 10,000 7000 
Wars arch bide Sich eine td 20.2 22.8 | 1.209 1.28 oe | 10.7 
yi’...| 0.65 | 0.62 |2.0X10°3/1.0XI10 4/6.3X10 5|1.9X10 °|6.4X10 4|2.3 X10 © 
Wie tr cace .05 .62 10.0 17.4 0.25 O.12 | 0.32 0.08 
AZ... <|"0:00:|0;00:'| =Fo.22 4-0-23' | =0:07 =6.11 —=0.05 —=©,14 


22. The maxima of spectral lines are complicated phenomena, 
often depending more on the properties of other elements than on 
those of the element that produces them. Moreover, they can rarely 
be determined with accuracy from existing observational data. 
Some better way of determining electron pressures in stellar atmos- 
pheres is to be desired. The most trustworthy data appear to be those 
for the general level of ionization in the solar atmosphere derived by 
the writer,” and in sun-spots, by Miss Moore.* ‘These depend upon 
a general calibration of hundreds of lines. The conclusion of Min- 
naert and Wanders” that the atmosphere above the spots is in 
radiative equilibrium fully justifies their use as independent data. 

26 Adams and Russell, Mt. Wilson Contr., No. 359; Astrophysical Journal, 68, 9, 1928. 

27 Mt. Wilson Contr., No. 383; Astrophysical Journal, 70, 54, 1929. 

28 Mt. Wilson Contr., No. 446; Astrophysical Journal, 75, 321, 1932. 

29 Zeitschrift fiir Astrophysik, 5, 297, 1932. 
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For the disk, 7 =5740° and J,=8.5+0.1. For this value of J, V.= 
N,, and y=2; whence K = 23, P=46+9. For the spots, T=4720°, 
f=7.0t0.1, P=28+ 7. 

We shall now assume, as the best convenient approximation, that 
the solar atmosphere contains in addition to completely ionized and 
neutral gases a single ‘active’ or partly ionized element. If the opac- 
ity is produced by ions, we have, by (83), 


uw’ Mogr,=C WT ~2G(A, T)P2u, . (120) 


Here yw’ is the mean atomic weight in chemical units and M,= 
1.649 X10 “4, the mass of this unit. For the sun, go = 2.74 X 104. Set- 
ting as before \= 4400, we find for the disk G(A, T) = 741, and 


p’t = 2.33 X10°55, P?uU,=0.050U2 . 


For the spot, G= 2400, yw'r =0.081n). 

So far, we have made no specific assumption regarding the com- 
position of the atmosphere, but we must do so to find w,. Since hy- 
drogen is in great excess, we may set €,=1 in (83), obtaining u,= 
5+f.,(a’y) —f.(y), so that w, is always between 3 and 1. About three- 
quarters of the metallic atoms in the sun* are of Mg, Fe, and Si, 
with a mean / of 7.8 volts. Most of the rest are of Ca, Al, Na, and K, 
and are almost completely ionized, even in spots.** We therefore 
adopt a’=0.25, /=7.8, and find for the disk, K =95, y=0.49, u2.= 
0.58, u'r =0.029; for the spots, K =1.96, y=14.5, u.=0.61, pT = 
0.050. Since hydrogen is in great excess, u’ must be nearly 1. Even 
so, the computed values of 7, are about one-tenth as great as might 
have been expected. No permissible changes in the assumed solar 
data can do much to relieve this discordance. This may arise from 
imperfections of the approximate theory ;** but it is just what was to 
be anticipated if neutral hydrogen atoms produce most of the opac- 
ity, and it may be regarded as evidence in favor of this hypothesis. 

i H. N. Russell, Mt. Wilson Contr., No. 383; Astrophysical Journal, 70, 56, 1929. 

°C, E. Moore, Mt. Wilson Contr., No. 446; Astrophysical Journal, 75, 322, 1932 
(Table VIII). 

32 Professor D. H. Menzel has very kindly informed me that his unpublished work 


indicates that the value of 7 given by Pannekoek’s formula should be considerably 


increased. 
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Certain additional tests are available. Miss Moore, comparing the 
disk and spot, finds, in our present notation, 


>/ yl 


P =0.60+0.10 , V =1.70+0.11, 


the primes referring to the spot. Now, by (80), 
u’M.gN=aPu. (121) 


With the data already adopted, w’ = 3.36, w=1.33, whence P’/P = 
0.40.V’/N =0.68. The agreement is within the probable error.*3 For 
the number of atoms per square centimeter above the photosphere 
we find uw’ N = 1.02 X 10°". This does not include the non-ionized ma- 
terial. From direct calibration of solar lines, with the aid of Unsdéld’s 
results, the writer finds 8 x 10”. 

These results depend only on the assumption of an isothermal at- 
mosphere, and not on the law of opacity. The latter determines the 
change of P and N with temperature. If u=constant, 7,=constant, 


oe eT? 
Pe GQ, T)’ 


If r,=constant, 


vt: (123) 


where a’(1—e) =e’, the fraction of the atoms which are completely 
ionized to the lowest temperatures considered. 

With the values just derived for the sun as standards, those for 
any other star may be computed. Since w is a function of y, and 
therefore of P, the solution must be made by approximations, which 
converge rapidly. 

On introducing the solar values, (122) gives 


log P=+0.16+3 log : +7 log r—} log u.—} log G(A, T) . (124) 


33 Cf. the more detailed calculation (loc. cit.). 


34 Mt. Wilson Conitr., No. 383; Astrophysical Journal, 70, 63, 1929. 
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Then, for the spot, 


log P’=1.23+4 log u;. 


If we assume for the first approximation that #;=0.58, as for the 
disk, log P’ =1.35, P’ = 22.2; but K’ =1.96, y’=11.4. Hence, «= 
0.616, which gives P’ = 21.6, y’=11.1. This yields ui =0.618, P’ = 
- P, b] 
21.6. The third approximation is superfluous. 
We now have w’ = 3.20, and by (121), V’/N=1.12. With atomic 
opacity, 


PV u,« T (eg)? ; (125) 


For the disk, #,=0.179. For the spot, y is large; and, as a first ap- 
proximation, we may set u;= 3; whence, P’=20.5, u,=0.431, 
P’=22.0. The third approximation gives P’= 21.8, yW=11.1, u4= 
0.442, and N’/N =1.14. 

The agreement of the results of the two assumptions is a coinci- 
dence. The computed values of P agree tolerably with those ob- 
served; those of V’/N are seriously discordant. With a large value 
of a’, that is, with more of the easily ionized metals, the agreement 
on the second hypothesis would be better. 

23. We may now determine a’ from the observed maxima of Fe+" 
and 7i+ (Table XII). For various assumed temperatures of maxi- 
mum, we calculate P! uw, from (124), taking log (g/g,.) from (118) 
for giant stars, and also K for hydrogen. Then from (119) we find the 
values of a’y which make the line under consideration a maximum at 
this temperature. The values of ., y, and a follow at once. If the 
opacity is due to neutral atoms, we may use (125), remembering 
that, on our assumptions, e’ =a’ when the hydrogen is being ionized, 
but e’ =}a’ inthe sun. We then obtain the results given in the lower 
part of Table XII. 

The observed maximum for Fe* is in class Fs, for which, by Table 
VIII, s=1.88; whence, if the opacity is due to ions, a’=0.001. For 
Ti* it is at F8 (s=2.05) and a’=0.0005. With opacity due to neu- 
tral atoms, a’ =0.0004 in both cases. 

A change of 0.2 spectral class in the assumed position of maximum 
(which is by no means improbable) changes a’ by a factor of 4. 
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There is also some uncertainty regarding the value of gravity, but 
this is less serious. According to Strémberg, super-giant stars of 
classes F and G average about four magnitudes brighter than nor- 
mal giants. If they obey the mass-luminosity law they should have 
about six times the mass of normal giants and eight times the di- 
ameter (allowing something for lower temperature), so that g should 
be about one-tenth the assumed value. If all the observed stars 
were super-giants, the tabular values of a’ would have to be multi- 
plied by about 3. Any permissible correction must be much smaller. 

The conclusion that a’ is of the order of 0.001, that is, that hydro- 
gen is about a thousand times as abundant as all the metals together, 
appears to be well founded. 

The strength of the Balmer lines in the sun indicates a still greater 
abundance of hydrogen, but, as is well known, the determination is 
vitiated by the Stark effect. Unséld, by comparison of the intensi- 
ties of the (K) line and the continuous absorption beyond the Bal- 
mer limit (which should be free from disturbance by the Stark ef- 
fect), finds for class A2, log, Viz/ Neu = 4.2, but reduces this to 3.0 by 
an empirical allowance for departure from thermodynamic equilib- 
rium. Since calcium contributes about 4 per cent of all the metallic 
atoms,*> this makes a’ =0.0016 without the empirical correction, and 
0.025 with it. The reality of the “correction,” which was suggested 
four years ago by the writer, now appears very doubtful. 

Menzel’s careful study of the chromosphere* leads to a’ = 2X 10~°. 
This value involves an empirical correction for departure from ther- 
modynamic equilibrium. 

24. The relative absorbing power of ionized and neutral hydrogen 
atoms may now be estimated from the maxima of H or Mg* by (106) 
and (111) with €,=o: 


a’ug(ti tr.) =C3T > 2 PP; u.,=Hu,+ An; . 


Now bi =a'A. If bya’ and a’y is finite, we may set u,.,=;. But 
= 3 3 
if y is finite, 
pe aaah u;~a'(1+3y) , 


5 H. N. Russell, Mt. Wilson Contr., No. 383; Astrophysical Journal, 70, 56, 1920. 


36 Lick Observatory Publications, 17, 281, 1931. 
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so that 


We then have 


where b’=b/(1+b). Since b’y is determined by the condition for the 
‘ oe a eo 
maximum, the last expression is approximately known. If a’ has 
been determined as above, y may be found for any given tempera- 
ture, and then 0’, then a better value of y, and so on, with rapid con- 


vergence. 


Table XIII gives the results for giant stars, on the assumption 
that a’=0.0004. When y is large, b’<a'. Computation shows that 


this relation is very nearly true up to a’=0.04. The uncertainty of 
the precise spectral type at which the maximum occurs, and that of 
the actual temperature of A stars, is too great to permit even an 
attempt to evaluate 0; but it is clear that there is abundant latitude 


Pu; = a’ K?(y?+ 23) a Tg 
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un; _ I+b+3 by 
b+ 5 by © 


U3 


TABLE XIII 


iH 
y b 
4.63 10.21 
15.1 054 
| 05.7 O.OI1! 


to satisfy our assumed condition, a’«b<«1. 


25. The degree to which existing theories agree with observation 
cannot be well tested by computations from assumed maxima. It is 
much better to start with some definite assumption simple enough 
to be workable, compute the electron pressure ‘‘at the photosphere”’ 
and the quantity of gas above it by the formulae already given, and 
then calculate the number of atoms which should be at work in a 
given line. Since the importance of opacity due to neutral hydrogen 
atoms is still uncertain, it is assumed that the opacity arises entirely 


la 
1+d'y ’ 


Met 


i 


5-15 
10.9 
74.4 


b 


0.38 
oSS 
0.018 
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from encounters of ions and electrons, and that, at a given tempera- 
ture, only one element is “‘active’’ or partly ionized, the rest being 
completely ionized or neutral. Our fundamental equations then are 


) 


i 
(1—€)ugv =P | Thin pm; (126) 


log PV u,=0.16+4 log : +?log7—}logG(a,T), (124) 


u,=3(1+a' —a’e,) +f.(a'y) — (a +46) f2(y) , (127) 

log 8 =0.76—0.973 (giants) , (1176) 
20 

= —0.65+ 0.28 (main sequence) , (117@) 


z=11,600/7; y=P/K. Logarithms are here to the base 10. Equa- 
tion (So) has been modified so that V represents the number of all 
atoms, not that of atoms of the active element, as before. €, is the 
fraction of non-ionized atoms in the whole atmosphere; a and a’ 
those of partly and completely ionized elements among the rest. 
Equations (117a) and (1176) are empirical ($18); the numerical con- 
stant in (124) depends on solar data ($22). When most of the ions 
come from the hydrogen, G (A,7) in equation (11) should be replaced 
by H(T), equation (12). The composition of the atmosphere is as- 
sumed to be 999 parts hydrogen to 1 of metals (a’=0.001). Of the 
metals, three-fourths are iron, three-fourths of the rest sodium, and 
the remainder potassium. This gives a tolerable representation of the 
composition of the sun, without being too complicated to work with. 
Helium is neglected, since, according to Unséld’s determination (the 
best yet available) ,‘’ its abundance is only 1 per cent that of hydro- 
gen. Oxygen has so nearly the same I.P. as hydrogen that it may be 
combined with it. Second ionization is neglected, since it affects the 
whole number of free electrons by 1 per cent at most. The field of 
temperature falls into several ranges, in each of which one of the ele- 
ments undergoes ionization. Within each range, a’ and e¢, are con- 
stant, so that « and u, are functions only of y, and they or their 


37 Zeitschrift fiir Astrophysik, 3, 99, 1931. 
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logarithms may be plotted, and the values of A (for the active ele- 
ment) and G (or //) tabulated for various values of 7 or z, as also 
that of log P|! u.. We then have log yl #, at once, and log y almost 
by inspection, and then log P and log NV. 

The results for the main sequence are shown in Figure 1, in which 
log NV is plotted against 3’=2 Mod = 5040/7. The curves marked //, 
Fe, Na, and K show how log NV or log P would vary in an atmosphere 
containing the assumed proportions of the given element in which 
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the other elements were completely neutral or wholly ionized under 
all conditions. The nearly straight parts of the curves, sloping up- 
ward to the right, correspond to temperatures at which the active 
element is practically all neutral or all ionized; the intermediate part, 
with the inflection, corresponds to its ionization. The final part of 
one curve and the initial part of the next are asymptotic to curves 
representing an atmosphere containing only completely ionized and 
completely neutral atoms in the appropriate proportions. These 
curves are slightly concave upward. The curves for H and Fe, and 
again for Fe and Na, approach very closely, showing that the ioni- 
zation of one is almost complete before that of the next begins. 
This is not the case for Na and K. The course of the curve when 
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both are undergoing ionization has been calculated with the aid of 
equations (54), (58). The other curves approach so closely that the 
transition can be drawn in freehand. The rise at the lowest temper- 
atures (s’ >1.8) results from the assumption that no other source of 
opacity is present besides the interaction of the few electrons and 
ions which come from the potassium. When 3’ = 2, only 1 in 105 of 
all the atoms is ionized, and opacity due to other sources probably 
exceeds that here considered. This would flatten out the rise at 
lower temperatures and might reverse it. 

A complication arises in the curve for hydrogen. With very small 
ionization almost all the free ions are metallic, and the function G 
must be used in the equation. With high ionization, the hydrogen 
predominates and // is to be used. The number of metallic ions is 
a’ N, while that of the hydrogen ions is .V,. The contributions of the 
two to the opacity are proportional to Ga’N and IIN,, or, by (80) 
and (81), to aa’Gu and Hu,. In the region of transition, where the 
two are comparable, G and H/ are not very different, and a linear 
interpolation between the values computed with each separately 
gives accurate enough results. 

The constant 2 has been added to log P in Figure 1 to get it into 
the picture. Here the curves for Fe and Va show shallow depres- 
sions due to the change of «, from o.50 for both high and low 
ionization to a maximum of 0.66, when y= 3.5. Those for AK and H 
drop lower, since uw, rises to 1. The net effect of all these influences is 
that log P varies almost linearly with z’ except at the highest tem- 
peratures. 

Figure 2 gives similar data for giant stars having the assumed re- 
lation between surface temperature and gravity. The change of log P 
with temperature is here very rapid, and still more nearly linear. 
The empirical expression log P = 3.46 — 2.742’ represents the values for 
3’=2.40 to 3’=0.50 with an average residual of only +0.025. At 
this temperature the giants run into the main sequence, and to use 
the formula beyond this limit is illusory. The values of V for the 
giants show an even greater drop, while the hydrogen is being 
ionized, and increase gradually toward lower temperatures, while 
those for the main sequence decrease. It is now evident why the 
assumption of an atmosphere of constant extent gives in many in- 
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stances a tolerable approximation. The mutual interaction of op- 
posing factors causes the amount of material “above the photo- 
sphere” to vary relatively little from zs’ =0.8 (1 =6300°) to the low- 
est temperatures. The amount is again nearly constant from 10,000° 
upward. The rudimentary theory based on a constant opacity co- 
efficient and extent of atmosphere therefore gives a fairly good repre- 
sentation of the effects of temperature on line intensity from classes 
M to F5 and again from Bog upward. It fails in the intermediate re- 
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gion and does not account for the absolute-magnitude effects even in 
the cooler stars. 

26. We are now in a position to calculate the number of atoms 
engaged in the production of any given line. The number of atoms of 
the given element per unit area above the photosphere is NV” = Ne”. 
For those of neutral or ionized atoms, we may use Milne’s approxi- 
mation and set Nj’=N”y/(2+y), Ni’ =2N"/(2+y). The error of 
this approximation is unimportant for our present purpose, except 
when ionization is just beginning in the presence of other fully 
ionized material, that is, for ultimate enhanced lines near their 
first appearance, and no such lines are observable, unless perhaps 
in the fainter red dwarfs. When second ionization is sensible, 
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we have N=N,(1+3y,1+3yy.). Since y,>y., we may set N= 
N.(1+391)(1 +49.) with no sensible error. 

The number of atoms in the appropriate excited state is M¢’= 
Niqe-"*/B; but only a certain fraction, f, of these will be engaged 
in the production of a given line. For the leading members of strong 
series, f is nearly unity; for improbable, though permitted, transi- 
tions, it may be very small. Values have been calculated by wave 
mechanics in a few of the simplest cases, such as the Balmer series.* 
For other strong lines it may often be estimated, e.g., by Kronig’s 
rules for intensities in triads. For weak lines it may be found em- 
pirically. If Z is the number of atoms actually at work on the line, 
we have for arc lines 
Vo 


—Ez’. 
2+ y. 


logy Li’ =log N+log e’’ +log W 4 tog a +log 
B 2+%: 

The abundance factors €’’, which have been taken from the writer’s 
work on the sun’s atmosphere, are proportional to the quantity called 
T in Table XIV.*° From this table, log 7 =8.1 for all the metals (in- 
cluding silicon). With a’=o.001, log 7 =11.1 for hydrogen. Then 
log e=log 7—11.1. The partition functions B have been taken as 
constant, which again introduces no significant error. The weights q 
apply to the individual lower energy state for the line, except for un- 
resolved blends, when the whole weight of the blend is taken. Reasons 
for the assumed values of f are given in the notes to Table XIV, 
which summarizes the data for a number of important lines. 

Figure 3 shows the predicted behavior of log Z for the metals and 
hydrogen in the main sequence, and Figure 4 in giants. The shapes 
of the curves reveal the processes which are at work. For example, 
the curve for Fet in Figure 4 rises steeply at low temperatures (near 
z’=1.3) because the iron is rapidly becoming ionized. The smaller 
slope beyond z’ = 1.1 corresponds to the Boltzmann factor, practically 
all the iron being ionized. The maximum at o.8 and the ensuing drop 
result from the increase of opacity due to ionization of the hydrogen. 
As this comes to an end, the curve eases off and would turn upward 
again were it not for the second ionization of the metal, which causes 
38 Y. Sugiura, Scientific Papers, Tokyo, 11, 1, 1929. 

39 Mt. Wilson Contr., No. 383; Astrophysical Journal, 70, 56, 1929. 
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the final rapid drop. The predicted behavior of the metallic lines 
agrees closely with observation. As an illustration, the maxima 


TABLE XIV 


DATA FOR INDIVIDUAL LINES 


El d log € I Bo B, I q log is Notes 
(qf/ B) 
I] 4340 0.0 | 13.54 2 I 10.15 8 |—1.35 |—0.75 I 
He 4471 |—2.0 | 24.48 2 2 20.87 8 |—1.2 |—0.6 3 
Tle 4686.) 0.3... 54.16 2 I 43.16) I6 |=0..25 |-FO:7: 2 
O 7792 ('—=2.02| £3150 9 4 Q.1! 3 {0.1 -0.3 2 
O MOA sic 5 acs 34.904 4 9 22.90 6 |—0.6 |—0.3 5 
Ug 5183 |—-3.3 7.69 2 2 2.67 5 |—0.6 |—o0.2 4 
Vg AAST |. 06s 14.97 2 I 8.83} 10 |—o.1 |+0.6 2 
Si .s| 3905 |—3-8 8.12 fa) 6 1.90 I —O.15:/—I.1: 2 
oe SENDS asec eats 16.27 6 2 9.79 6 |=—0.2 -0.2 2 
Si ee a Eee 33.30 2 2 18.92 3 —oO.! 0.0 2 
S1 4OSS |...... 44.95 2 I 23.95 2 |=0.1 =5, 1 2 
Ca 4226 |—4.4 6.09 2 2 0.00 2. \=O.2 =O). 2 
RM teenies siete tee aie or 1.385 5 —fi0 [=r 3 
Ca* Bee Po Eee 2 I 0.00 2 |—0.2 |—0.2 2 
Fe. 4383 |—-3.9 7.83] 45 60 1.48 9 |-0.8 |—1.5 6 
RF 5 is ARR ale ne ries 16.9 60 30 2.80} IO |—2.3 |—3.1 7 
Zn ..| 4810 |—6.2 9.36 2 2 4.06 5 |—0.6 -O.2 4 
SAS eee Penareee) fame sa) TI28O 2 Doom (resteaeteds belat carted rachis ais rs sree (SOA 


NOTES TO TABLE XIV 
t. Calculated from wave mechanics. 
2. First line of strongest series from the given limit. 
3. Second line of strongest series. 
4. First line of weaker of the two subordinate series. 
5. Leading line of strongest multiplet of principal triad. 
6. One of four strong lines of nearly equal intensity from this level. 
7. Combination between terms having different limits in Fett. The strong multi- 
plets from this level are in the far ultra-violet. The total intensity of the observable 
solar lines from this level corresponds to log = 2.6 as against 4.6 predicted for all lines. 
Hence log f is about —2.0 for all the visible lines and — 2.3 for the strongest from the 
given component level. 

For the higher ionizations, the columns headed By, and B,; give the values for the 
lower and higher states involved. 

In this table the values of g for states involving two equivalent s-electrons (such as 
the normal states of He, Mg, Ca) have been doubled. I am indebted to Professor E. U. 
Condon for information that recent developments in wave mechanics show that this 
should not have been done. As it was done, however, in the calculations of the abund- 
ance of the elements in the sun, on which the present computations are based, the data 
are at least self-consistent and satisfactory for illustrative purposes, which is the ob- 
ject of the present diagrams. A recalculation would affect the computed maxima by 
much less than the errors of observation. 
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for the giant series are shown in Table XV. The values of 2’ have 
been translated into spectral type with the aid of Table VIII. 
The observed values, from Miss Payne’s Stellar Atmospheres,” 
usually agree as well as might be anticipated from the character 
of the data, except that Hy and Mg* come at too low a tempera- 
ture (see below). For Cat, on the other hand, the maximum in 
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FIG. 3 


class K, which has hitherto been a great puzzle," is explained. 

What was overlooked in the earlier theories was simply the fact that 

the maximum for an ultimate enhanced line for an atmosphere of 

uniform depth is very flat-topped. For the giants, the diminished 

pressure and opacity at low temperatures throw the maximum to the 

point where first ionization is nearly ended; for the main sequence, 
0 Pp. 127-130. 


4« Cf. Milne, op. cit., 92, 183, 1932. 
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the opposite change in gravity throws it to the onset of second ioniza- 
tion, in class F5, with very little change down to Ko. 
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The relative values of log L for z =0.88 are of course in agreement 
with the solar data from which they were derived; but at other tem- 
peratures there is unforced agreement with observation; for example, 
the approximate equality of the maximum values of log L for Fe* 
and Mg*. Miss Payne finds a maximum intensity of 6.0 for the lat- 
ter and of 9.0 for other Fe+ lines comparable in strength with the 
former. The latter figure must be raised by the inclusion of super- 
giant stars. 
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For hydrogen the well-known and enormous discrepancy appears, 
the computed maximum L for //y being only ten times greater than 
for \ 4481, while the actual strength indicates a value fully a thou- 
sand times greater. This can hardly be attributed, even in part, to 
an abnormality in the Boltzmann factor, since \ 4481, for which E is 
seven-eighths as great, agrees very well with prediction. It is doubt- 
less due to the “pressure effect” of the fields of neighboring ions, as 
discussed in detail by Unsold. The values for calcium in giants and 
the main sequence are compared in Figure 5. The enhanced line (K) 
is always stronger in giants. At low temperatures, the difference 
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should be very conspicuous, and it is actually so in class M. The arc 
lines also are predicted to be stronger in giants than dwarfs, below 
3’ = 1.2, corresponding to gG8 and dK3. For temperatures above this 
value the intensity should be greater in the dwarfs, if the comparison 
is made between stars of the same temperature; but for stars of the 
same spectral class the temperature difference very nearly counter- 
acts this inequality. This prediction could well be tested by observa- 
tion with high dispersion. It has, of course, no bearing on determina- 
tions of absolute magnitude from the relative intensities of neighbor- 
ing lines of different sorts. 

The predicted maxima for Fe+ and Mg* are much too near to- 
gether. It was shown in $20 that the assumption of opacity due to 
neutral atoms would correct this. By (106) the total opacity is thus 
proportional to Gu,+ Au, (when metallic ions predominate). In $22 
it was found that the term Gu, accounted only for an opacity ut = 
0.029 in the sun, while the actual value is probably about 0.4. With 
G=741, u,=0.58, u;=0.18, this involves A=31,000. Now A= 
C./C,é and, with our assumed data, e’= 1/4000. Hence C./C,=7.9. 
The opacity due to hydrogen ions alone is proportional to 80 C,; that 
due to atoms, to C., so that to represent the solar data, we must 
assume the latter to be one-tenth of the former. This looks like too 
high a value, but it is at least interesting as a typical case for calcula- 
tion. We find, instead of (124), 


log P=0.48+3 log 8 +? log T—} log (Gu.+ Au;) . 
So 


The calculation of P and NV then proceeds as before, with very little 
more labor. The values of log P and log N for the main sequence are 
shown in Figure 6. .V is now very nearly constant for the cooler stars, 
but diminishes considerably less for the hotter ones than before, 
while P increases rapidly during the ionization of the hydrogen. The 
increase in opacity, which on the former assumptions was due to 
both the electrons and the ions, is now thrown on the electrons alone, 
which therefore have to be more numerous. 

The values of log L for Hy, Mg*,and Fe* are shown in Figure 7, 
on both the new assumption and the old. The maximum for Fe* is 
hardly shifted at all, while those for Hy and Mg? are shifted by 
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—o.12and —o.18, or by 0.22 and 0.33 spectral classes. This extreme 
assumption removes less than half the observed discordance. The 
rest may be attributed to distortion of the energy-curve. Curves for 
Si,O, and He in various stages of ionization for main-sequence stars 
are shown in Figure 8. The double maximum for Si* is produced by 
the ionization of hydrogen. 
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In Table XVI the computed maxima are compared with Miss 
Payne’s observations. The positions of the maxima agree closely with 
observation, which is not new. The predicted maximum intensities 
for Si*, Si**, Si***, and O* are nearly equal, as are the observed 
values. The Si line at \ 3905 is not so strong in proportion as 
might be expected; but the value of f for this line is uncertain. The 
results for O (\ 7772) are rough, since the values of G and // are 
different at this wave-length. The solar intensity of this line is 2; of 
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d 4810 Zn, 3; and of \ 4481 Mg*,o. Since the values of log L for the 
three lines are —5.5, —5.2, and —5.8, the agreement is tolerable. 

The abnormal strength of the lines of H and He appears to be con- 
fined to these elements. Measures of line contours and effective 
widths in spectra of classes A and B would be of great value in test- 
ing the theory. 

27. The foregoing discussion leaves the maxima of stellar spectral 
lines in a satisfactory state. The assumption that the distortion of 


_ the energy-curve, which is known to exist, makes the color-tempera- 


tures of A stars too high removes the last serious discrepancy. 


TABLE XVI 


| 

















El. 3! te | Pred. Sp. | Obs. Sp. | Obs. Int. log L+10 

Si. 20434 5300° G2 | Go mee 7.9 
Hy 67 7500 | A7 | Ao 16.0 6.7 
O.. .68 7400 A8 | EN eee NS 5-9 
ore f 68 | 7400 A8 | ice 

St. | 43 | 11700 BS} Ao 9.3 $e 
He. . at 16200 | B4 B3 11.8 Ce 
RY uaa 96 20000 | B3 Bi 10.0 ee 
OF =. .24 21000 | B2 Br 12.0 4.7 
Si .19 26000 | Bo Bo Q.2 4-5 
TRG Ys. 0-05 0.14 36000 | O7: O >6 4.2 


Two grave discordances between theory and observation remain. 
One, the abnormal strength of the lines of H and He, may be attrib- 
uted with some confidence to the Stark effect, which is probably ca- 
pable of accounting for it quantitatively. The other, and really the 
greater discordance, occurs in the spectra of red giants, where the 
Balmer series, and the enhanced metallic lines as well, are enormous- 
ly stronger than the present theory predicts. The Boltzmann factor 
has little influence in the latter case, and it seems to be impossible 
to account for the facts on the assumption of an isothermal atmos- 
phere. 

If, however, the opacity in the wave-length considered were much 
smaller than the Rosseland mean, we would receive light from deeper 
and hotter layers. An exact study of ionization and excitation in 


42 Adams and Russell, Mt. Wilson Contr., No. 359; Astrophysical Journal, 68, 27-35, 
1928. 
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such an atmosphere, in which the radiant flux is not of black-body 
type, would be difficult. It is easy to see, however, that the number 
of highly excited atoms would be relatively increased. This effect 
may account for the curvature of the plots of log M against E found 
by Adams and the writer,* thus making the suspected departure 
from thermodynamic equilibrium a simple difference in temperature. 
Part of this effect, however, may be due to inaccuracies in the cali- 
bration of the very strong lines of low £ in the red stars; but this cir- 
cumstance cannot account for the behavior of the enhanced lines. 
Whether such a difference between the monochromatic and Rosse- 
land opacities exists demands theoretical investigation. One would 
not expect to find it in stars where the maximum of the energy-curve 
occurs in the observed region, and there is little or no evidence of it 
there—the diagram given by Adams and the writer“ being repre- 
sented as well by straight lines as by curves. It is strongest in Betel- 
geuse and Antares, where the energy maximum is well in the infra- 
red. 

The increased pressure in the deeper layers inhibits ionization, 
while the high temperature favors it. A given percentage change in 
pressure produces the same effect at all temperatures, but a change 
in temperature is far more effective at low temperatures. The ap- 
pearance of this effect in class M is therefore intelligible. 

Present theory, however, indicates a high, rather than a low, 
opacity in the violet for such stars, and some other explanation may 
be preferable; for example, Unsdld’s* suggestion that radiation in the 
far ultra-violet beyond the series limits of the metals may come 
through strongly enough to produce the excitation and ionization. 
A complete theory of stellar absorption lines will clearly be a com- 
plicated affair. 

One more point may be mentioned in closing. Suppose that some 
stars contain a considerably less overwhelming excess of hydrogen 
than the average. If the difference extends to the interior, as well as 
the surface, these stars will be brighter than the mass-luminosity re- 
lation indicates; they will be, or tend to be, super-giants. The maxi- 
mum of the Balmer lines will come at a higher temperature than 

4 [bid., p. 20. 


4 Ibid., p. 23, Fig. 2. 45 Op. cit., 1, 1, 1930. 
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usual, but, owing to the low density, the lines will be sharper and 
may appear fainter than in normal stars. The lines of other elements 
will be stronger than usual, compared with hydrogen, especially the 
enhanced lines of the metals. The arc lines may be fainter, owing to 
the high ionization. Lines of high excitation, ordinarily absent, may 
appear. 

These predictions of theory amount almost to a description of the 
spectrum of Alpha Cygni, and the suggestion that this and perhaps 
other c-stars are deficient in hydrogen appears plausible. To attrib- 
ute the deficiency to partial exhaustion of hydrogen by processes of 
atomic synthesis in these very luminous stars is tempting. 

PRINCETON, N.J. 
May 2, 1933 














NOTES 


REMARKS ON DARK NEBULAE 
ABSTRACT 


An attempt is made to explain the shapes of dark nebulae. It is suggested that the 
larger particles of a nebula will concentrate near the stars, while the smaller particles 
will be repelled into more distant regions. 

We have no definite indications of the size of particles forming the 
dark nebulae of Barnard and of Wolf. The criterion of mass required 
to produce the extinction of light observed in dark nebulae rules out 
gaseous matter as their chief constituent.’ Debye’s theory of the 
scattering of light by small particles (especially for non-transparent 
dielectrics) allows us to find this mass, given the extinction, the 
dimensions of the nebula, and the assumed size of the particles. I 
have performed such calculations, using Pannekoek’s data on the 
obscuration in Taurus and allowing 10° as a rather conservative 
upper limit for the total mass of this cloud. The results are not 
definite. The cloud in question may be composed of particles of all 
sizes, from 7.1077 to 10-* cm in diameter, capable of producing all 
kinds of scattering (A~4, At, A°-laws); a minimum mass of approxi- 
mately 15 corresponds to diameters of about 10-5 cm. The existing 
data on the color effects in dark nebulae are meager and somewhat 
ambiguous; further observations are greatly needed in order to de- 
cide the question. Nevertheless, even at present we can derive some 
conclusions as to the nature of dark nebulae by considering the 
forces acting on them. 

In a recent paper? the writer has shown that on account of radia- 
tion pressure small ideally reflecting spheres producing scattering 
according to \4or A‘ cannot be retained as permanent satellites in 

' Pannekoek, Proceedings of the Amsterdam Academy, 23, 720, 1920; Eddington 
Internal Constitution of the Stars, p. 388. 

2 Zeitschrift fiir Astrophysik, 4, 265, 1932. 
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the vicinity of supergiant or giant stars. For stars of classes Bo-Bs5, 
the lower limit of the radii of such particles turns out to be of the 
order of 0.2-0.04 cm, while late-type giants can retain particles of 
1o%-ro 4 cm in diameter. On the other hand, dwarf stars exerting 
much smaller repulsive forces can possess, in the form of satellites, 
very small particles capable of producing selective scattering.’ As 
to the order of magnitude, these results are also valid for non- 
transparent dielectric particles as well as for absorbing spheres of 
appreciable refractivity. 

The different behavior of small particles in the vicinity of giant 
and dwarf stars seriously affects their distribution in galactic space. 
General as well as local fields of radiation and gravitation cause a 
special distribution of dark matter which could be calculated if these 
fields were known. 

The properties of the general galactic field gradient represent the 
most difficult problem. At the present time we can only make some 
rough and indefinite guesses as to its direction and size. The accumu- 
lation of stars of high luminosity in and near the galactic plane gives 
some reason to think that there exists an appreciable general field 
gradient perpendicular to the galactic plane, tending, in higher 
latitudes, to drive away small particles. We have no knowledge con- 
cerning the projection of this gradient on the galactic plane, except 
assumptions by analogy based on the occurrence of dark belts in the 
equatorial planes of spiral nebulae. In contrast to the general 
galactic field, the local fields are more definite and can sometimes 
be inferred from the observational data. 

The most outstanding examples of the repulsive action of stars 
of high luminosity on cosmic dust are the ‘‘cometary nebulae.” 
These are feebly luminous, nebulous objects of comet-like shape with 
a more or less conspicuous star at the apex. At the present time the 
following nebulae of this kind are known: 


N.G.C. 2261 (Hubble’s variable nebula): Nucleus—R Monocerotis (spectrum 
nova-like). Nebular spectrum continuous. 


3 These calculations have been repeated recently by Schénberg and Jung (Astro- 
nomische Nachrichten, 247, 413, 19033), who reached essentially the same conclusions. 
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N.G.C. 6729 (variable): Nucleus—R Coronae Australis (Gpe). Nebular spec- 
trum continuous. 

N.G.C. 2245: Nuclear spectrum Br. Nebular spectrum continuous. 

N.G.C. 1555 (Hind’s variable nebula): Nucleus—T Tauri (spectrum peculiar). 

Anonyma (a= 63™, 6=+18°42"; 1920): Discovered by Hubble. Nuclear spec- 
trum Br. Nebular spectrum continuous. 

Anonyma: Nucleus CPD — 25°2286. Spectrum Ma. Discovered by Van den Bos. 

AAnonyma: Nucleus B.D.+ 28°645 (RY Tauri, irregular variable, F8). 


From a purely mechanical point of view, these cometary forma- 
tions may be explained in much the same way as ordinary comets. 
A dust cloud containing particles of various sizes, small as well as 
large, if approached by a very luminous star, necessarily takes the 
shape of a comet, small particles being repelled, while the envelope 
of their orbits assumes the shape of a paraboloid of revolution. On 
the other hand, large particles (if present in sufficient quantity) are 
attracted by the star from the head, both head and tail shining with 
reflected starlight. Such a nebula is in general doomed to dissipa- 
tion. In some cases, when the head exists and is sufficiently opaque, 
the large particles which compose it may produce a screening effect, 
partially protecting the smaller particles. 

In four cases out of seven the nuclear stars are variables (two of 
them are nova-like). For purely physical reasons it is difficult to 
regard the variability of the nuclear stars as being caused by the 
clouds. On the other hand, nova-like outbursts would produce more 
or less symmetrical gaseous nebulosities. The phenomenon under 
consideration can be better explained in the following way: Stars 
showing nova-like spectra, and some irregular variables, are prob- 
ably the most luminous objects in the sky. Brighter supergiants 
produce sharper cometary forms which are better illuminated and 
therefore easier to discover. The nuclear star, if seen through the 
axis of the nebula, would have a nebulous appearance. If the star 
is behind the nebula and the nebular tail, a color excess is to be 
expected and, as we know, is actually observed in many nebulous 
stars. If the nebula is optically thick and is seen edgewise, the 
nuclear star would be invisible and the nebula itself would be pro- 
jected on the background stars as a dark, comet-like marking. Such. 
probably, are Barnard’s ‘‘dark comets,” B64, etc. 
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We may present the problem in a somewhat different manner. 
Consider a group of early (or late supergiant) stars initially im- 
bedded in a dense dust cloud composed of large as well as small 
particles. What will be the ultimate fate of such a cloud? We can 
expect that the larger particles will concentrate around individual 
stars as their Keplerian satellites, while the small particles will be 
dissipated, forming some kind of outer structureless obscuration. 
Such a steady state is probably attained in the Pleiades. On the 
other hand, under the protection of the large particles, if they form 
optically thick layers around the stars, small particles can be re- 
tained in a quasi-static state in the external parts of the nebula. If 
the steady state is not reached, these small particles can produce 
dark lanes and patches “radiating” from the central bright nebula 
which envelops the stars. As an example we may mention the dark 
and bright nebulae involving a group of early type stars near p 
Ophiuchi. These stars (p Ophiuchi, 22 Scorpii, B.D.—24°12684) are 
involved in bright nebulosities, showing continuous spectra with 
long dark lanes radiating from them in various selected directions. 

If concentrated bright nebulosities enveloping the stars consist of 
large particles, as is here suggested, they cannot give anomalous color 
effects in their reflected light. Recent observations of relative color 
indices, star-nebula, did not reveal any differential color effect,‘ and 
thus corroborate these conclusions. 

Another example of the repulsive effects may be suspected in 
Br44. This splendid dark nebula has a parabolic shape; its axis, if 
prolonged on the celestial sphere, passes directly through y Cygni 
(and through a Cygni farther on). The truncated summit of the 
nebula is 4° off y Cygni; the nebula itself extends up to 4° in the op- 
posite direction, fading out and diffusing into the regular stellar back- 
ground. Assuming that for y Cygni (F8p), according to a determina- 
tion made at Victoria, M = — 3™, we find the distance of the nebular 
summit from the star to be about 7 parsecs (if the axis of the nebula 
is tangential to the celestial sphere). 

If the local field is not central but results from several dispersed 
stars, its effect on a dark nebula should be very complicated; in this 


4Q. Struve, C. T. Elvey, and P. C. Keenan, Astrophysical Journal, 77, 274, 1933. 
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case the equipotential surfaces can take any imaginable shape. If 
such a field is repulsive, its strength would be proportional to 
er *, r being the optical depth reckoned from the outer boundary. 
Such a field would result in a one-sided compression of the cloud, 
producing one sharp boundary, the opposite boundary being diffuse. 
Every student of dark nebulae notices a remarkable property of 
these objects which has been described by Barnard as follows: “In 
many cases one side of a dark marking is very definite, while the 
other side is diffused. This occurs so often that there must be some 
reason for it.’’> Striking examples are B50, 143, 160, etc. It may be 
concluded, therefore, that in some cases the outer boundaries of 
dark nebulae represent equipotential surfaces, while in other cases 
they can be interpreted as the envelopes of trajectories or lines of 
flow. 

It goes without saying that the action of repulsive forces can 
explain only part of the observed contours of dark nebulae. It is 
probable, however, that dark nebulae represent rather ephemeral 
formations whose existence and evolution are determined chiefly by 
the local field of radiation and gravitation produced by stars. 
Internal forces are of minor importance. 

B. P. GERASIMOVIC 
POULKOVO OBSERVATORY 
August 1933 


VARIABILITY OF BRIGHT Ha IN a CYGNI 


ABSTRACT 


The emission line of //a in Deneb is found to vary in intensity. 


Some time ago we found that bright Ha in the spectrum of Rigel 
(8 Orionis) varies in intensity. Since this line resembles in structure 
the same line in a Cygni, we decided to observe the latter with our 
grating spectrograph, which is used in connection with the 69-inch 
reflector of the Perkins Observatory. The plates indicate a large 
variation in the intensity of the bright line, which may possibly be 


5 [bid., 49, 1, 1919. 
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associated with a variation of the absorption lines. The period may 
be in the neighborhood of several months, although a short period 


is not impossible. 


Date pig ool | Qual. | Date nr Qual. 

1933 Mar. 9.46... Weak | p 1933 Aug. 22.15...| Weak f 
Apr. 29.38... Fair | p ji Aug. 24.06...) Weak g 
May 18.38... Strong \ ee a Aug. 26.08...) Fair g 
May 30.38... Fair 1 # Aug. 26.19...| Fair g 
June 15.37... Fair bg Aug. 28.10...| Fair g 
July 31.14.. Absent | g || Aug. 29.17...| Fair g 
Aug. 19.18... Weak | f || Sep. 2.06...| Strong g 
Aug. 20.06... Weak be & Sep. 7.05...| Fair f 
Aug. 21.06... Fair | g | Sep. 7.10...| Weak p 
Aug. 21.12...| Fair | g | Sep. 18.08...| Very strong | g 
Aug. 21.18...' Fair g | Sep. 28.06...| Strong p 
Aug. 22.06... Fair be Oct. 5.05...| Strong f 
Oct. 12.01...| Fair g 





| 


The emission line was practically absent on an excellent plate 
taken on July 31 and was very strong on an equally satisfactory 
plate taken on September 18. On the last plate the absorption line 
of Ha is split in two by a narrow emission line which is much 
weaker than the emission line on the red side of the absorption 
line. 

O. STRUVE 
F. E. Roacu 
October 13, 1933 
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REVIEWS 


Atom and Cosmos. By HANS REICHENBACH. Translated by Epwarp 
S.ALLEN. New York: Macmillan Co., 1933. Pp.300. Bound,$2.00. 


In this book, which is an enlargement of a series of radio lectures de- 
livered from Berlin some three years ago, the author sets himself the task 
of presenting to the layman an account of the important advances which 
physics has made in the twentieth century, and of the changes resulting 
therefrom in the basic concepts of scientific thought. 

In general, the material is presented as simply and as precisely as could 
be expected in a strictly non-mathematical work. The author treats in his 
first section the concepts of non-Euclidian space and relativity of motion, 
together with their astronomical significance. Following this he discusses 
radiation, its apparent dual nature as wave-train and stream of corpuscles, 
and the reconciliation of these two theories in the developments of the last 
ten years. In the next section, various pictures of matter as made up of 
molecules, atoms, or electrons and protons are given. The author speaks, 
too, of the early pictures by Rutherford and Bohr of the interior of the 
atom. The last chapter of this section, on ‘““The Wave Character of Mat- 
ter,” is a fitting introduction to that which follows, in which the author 
speaks in everyday terms of the concept of probability, of Heisenberg’s 
uncertainty principle, and of the difficulties concerning fundamental con- 
cepts and world-pictures which scientific thought of the present century 
has had to face. 

Whether the author’s prefatory remark that his presentation ‘‘does not 
presuppose any knowledge of the kind taught in schools’ may be taken 
literally is doubtful; however, the book does accomplish its aim, to give 
“the physicist’s way of thinking, and a general view of the results of his 
research.” 

PAUL RUDNICK 











